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Two-loop scalar self-energies in a general renormalizable theory 
at leading order in gauge couplings 

Stephen P. Martin 
Physics Department, Northern Illinois University, DeKalb IL 60115 USA 
and Fermi National Accelerator Laboratory, PO Box 500, Batavia IL 60510 

I present results for the two-loop self-energy functions for scalars in a general renormalizable field 
theory, using mass-independent renormalization schemes based on dimensional regularization and 
dimensional reduction. The results are given in terms of a minimal set of loop-integral basis func- 
tions, which are readily evaluated numerically by computers. This paper contains the contributions 
corresponding to the Feynman diagrams with zero or one vector propagator lines. These are the ones 
needed to obtain the pole masses of the neutral and charged Higgs scalar bosons in supersymmetry, 
neglecting only the purely electroweak parts at two-loop order. A subsequent paper will present the 
results for the remaining diagrams, which involve two or more vector lines. 
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I. INTRODUCTION 



The next generation of high-energy collider experi- 
ments, including the Large Hadron Collider (LHC) and a 
future electron-positron linear collider (LC), should dis- 
cover the mechanism of electroweak symmetry breaking. 
They also hold the promise of very accurate measure- 
ments of the relevant parameters of that mechanism. If 
supersymmetry is correct, then the lightest neutral Higgs 
scalar boson mass should be measurable to an accuracy 
of order 100 MeV at the LHC and even better at an LC 
[l|. The mass of the top quark, which is known to be 
strongly coupled to the electroweak symmetry breaking 
sector, can be obtained with an accuracy of order 1 GcV 
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at the LHC and 100 MeV at an LC 0. These will 
easily exceed the accuracy obtained by theoretical calcu- 
lations at one-loop order in the Standard Model and its 
extensions. Therefore, calculation at two-loop order, at 
least, will be necessary for detailed comparisons of can- 
didate models with experiment. In minimal supersym- 
metry, for example, the higher-order corrections to the 
effective potential and to the mass of the lightest neu- 
tral Higgs scalar boson are especially important for two 
reasons. First, effects proportional to the relatively large 
QCD coupling first appear at two-loop order. Second, 
the tree-level scalar potential has a nearly flat direction, 
so that loop corrections involving larger mass scales are 
relatively enhanced. 

A great deal of progress has already been made on 
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calculations of two-loop corrections to masses and self- 
energy functions, including specific results for the top 
quark and the electroweak gauge bosons in the 

Standard Model [Tjj , the Higgs scalar bosons in the Stan- 
dard Model 01 an d supersymmetry 0-[2(J, and meth- 
ods for evaluating two-loop self-energy integrals in gen- 
eral HH-IH. 

Calculations of two- loop corrections can be difficult. 
Since we do not yet know with certainty what the cor- 
rect model of physics at the TeV scale is, or even a com- 
plete list of its degrees of freedom, a prudently flexible 
strategy is called for. Accordingly, I propose to compute 
self-energy functions for particles in a completely general 
renormalizable theory. These "once-and-for-all" results 
can then be specialized to any perturbative model, as 
future developments warrant. 1 A similar approach was 
used for the effective potential in ref. I will use ex- 
actly the same conventions and notations in the present 
paper. 

I choose to use a mass- independent renormalization 
scheme based on dimensional continuation with mini- 
mal subtraction, in which physically measurable quan- 
tities (physical masses, cross-sections, decay rates, etc.) 
are always outputs, and the input parameters are renor- 
malized running masses and couplings. One can then 
conduct a simultaneous global fit to all observables. 
This is in contrast to on-shell schemes, which use some 
preferred observables as inputs. The use of a mass- 
independent scheme means that one does not need to 
commit in advance to a particular choice of input observ- 
ables. Another possible advantage of mass-independent 
schemes is that the running input parameters can be 
evolved to higher mass scales (for example, the Planck 
scale, the scale of unification of gauge couplings, the 
supersymmctry-breaking scale, etc.) using the renormal- 
ization group, to compare with the predictions of vari- 
ous candidate models. For non-supersymmetric models, 
the traditional MS scheme [6j4| is used, while for models 
based on softly broken supersymmetry one may use the 
DR scheme |65| , which is based on dimensional reduction 
|6fi| but eliminates the appearance of unphysical epsilon 
scalar masses in the formulas. 

Consider a general renormalizable theory containing 
vector bosons, fermions, and scalar bosons, the latter la- 
beled by indices i,j,k, . . .. The main objective of this 
paper is to compute the two-loop scalar self-energy 
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where is the external momentum, and the metric is 
either of signature ( — h++) or Euclidean. Note that s is 
taken to be real with an infinitesimal positive imaginary 
part to resolve the branch cuts. The self-energy function 
is gauge-dependent, but can be used to obtain a gauge- 
invariant |S3-[z3 P°l e squared mass, defined as the po- 
sition of the complex pole, with non-positive imaginary 
part, in the propagator obtained from the perturbative 
Taylor expansion of the self-energy function about the 
tree-level squared mass |7lj . Explicitly, for a scalar par- 
ticle with tree-level renormalized (running) squared mass 
to 2 , the two-loop pole squared mass 



s k = M k - iT k M k 
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can be obtained iteratively as follows. First, the one- 



loop approximation to the pole squared mass, s k 
appropriate 2 solution of 
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Then, defining 
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the two-loop approximation to the pole squared mass sY 
is the appropriate solution to 



(2) 



Det (m? - 4 2) )*« + 



0. 



(1.6) 



Repeated indices are not summed in the last four equa- 
tions. The prime in equation (|1.5|) indicates a derivative 
with respect to the argument s, which is then replaced 
by to 2 ,. Note that [Hfc]y is not a function of s. 

In this paper, the contributions from diagrams involv- 
ing zero or one vector lines are included; the rest will ap- 
pear in a subsequent paper. Note that these diagrams are 
already enough to compute the full two-loop pole masses 
of all of the neutral and charged Higgs scalar bosons in 
supersymmetric extensions of the Standard Model, ne- 
glecting only the purely electroweak contributions (those 
that involve g 4 , g 2 g' 2 or g' 4 ). The explicit results for the 
particular case of the Minimal Supersymmetric Standard 
Model will also be presented in a separate paper. 

The rest of this paper is organized as follows. Section 
^specifies the notations, conventions and methods used. 
Section ITTT1 reviews the self-energy functions at one loop 
order. Section llVl reports the results for two- loop scalar 



1 Of course, the correct model of physics near the TeV scale 
might be non-renormalizable, or non-perturbative, or not even 
described by a four-dimensional field theory. But supersymme- 
try, at least, provides grounds for optimism. 



2 In this discussion, the "appropriate" solutions are the ones that 
continuously go over to m 2 as the loop corrections to the self- 
energy are formally turned off by reducing the three- and four- 
particle couplings of the theory to zero. 
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self-energies in the MS and DR schemes, with vector 
bosons assumed to have a generic non-zero mass. In sec- 
tion IY1 the corresponding results for massless vectors are 
obtained. Section IVII contains an illustrative example, 
featuring several non-trivial consistency checks. 



II. NOTATIONS AND SETUP 

In the following, I consider a general renormalizablc 
field theory, consisting of real scalars Ri, two-component 
Weyl fermions ipi, and vector fields A£. After expand- 
ing around the vacuum expectation values of the scalars, 
the tree-level squared masses of these fields are taken 
to be diagonalized, and are denoted by mf, m 2 , and 
m 2 respectively. This does not mean that the masses 
of the fermions are necessarily diagonal, however. The 
fermion mass matrices M IJ and Mjj = (M IJ )* appear 
in the numerators of chirality-changing fermion propaga- 
tors. They can be non-zero when the fermions ipi and 
if)j have the same squared mass. In this squared-mass 
eigenstate basis, the interaction Lagrangian terms have 
the form: 

C s = -^X^RiRjRk - ^ kl R t R J R k R l , (2.1) 
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with repeated indices summed over, and metric signa- 
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C FV = gfA^a^j, 



ture ( — h++). The squared masses m 2 , m 2 , m 2 , the 
fermion masses M IJ , and the couplings X^ k , A u , y IJk , 
gaij ( gabij ^ g abl ^ and gf J are all renormalized running pa- 
rameters. Raising or lowering of fermion indices implies 
complex conjugation of the coupling, so: 



gf = {gfT 



(2.5) 



(I have omitted possible three-vector, four-vector, and 
ghost interactions, which will not play a role in this pa- 
per.) For more discussion, see ref. |63j |. 

The propagators for the scalars Ri are diagonal, and 
are given for four-momentum by 

1^-2- (2-6) 



The fermions have diagonal chirality-changing propaga- 
tors, 

-' k '-' T (2.7) 
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as well as chirality- violating propagators that are not, in 
general, diagonal: 

-iM IJ -iMu 
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In a general R^ gauge, the propagator for a vector boson 
with squared mass m 2 and four-momentum fc M is 
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k? + m 2 m " V k 2 + nr; 
where I have introduced the notation 

C x f(x) = [f(x) - /(£c)] /x, 



(2.9) 
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with ^ = for Landau gauge and ^ = 1 for Feynman 
gauge. Results below for loop-integral functions involv- 
ing massive vector bosons will be conveniently written 
in terms of the C x notation. For massless vector gauge 
fields, one has: 

lim[£J^)] = (l-0/'(0), (2.11) 

corresponding to a propagator 

ir]^/k 2 + i{£- l^k" /{k 2 ) 2 . (2.12) 

Infrared divergences associated with massless vector 
bosons are dealt with by first computing with a finite 
mass, and then taking the massless vector limit. 

In preparing the Lagrangian as above, there is a choice 
to be made regarding the vacuum expectation values of 
scalar fields. One procedure is to expand around the 
VEVs that minimize the tree-level potential. This would 
mean that there are no tree-level tadpole couplings. How- 
ever, in that case, there can be non-vanishing tadpole 
graphs at one-loop order and beyond, which must be in- 
cluded in the self-energy functions. A different choice is 
to instead expand around VEVs that minimize the full 
loop-corrected effective potential. This means that the 
tree-level potential, in general, will have tadpole coupling 
terms: 
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However, expanding around a minimum of the full loop- 
corrected effective potential means that the sum of all 
tadpole contributions always vanishes, through whatever 
loop order we are working. Therefore, as long as we ex- 
pand around a minimum of the full effective potential, we 
can safely neglect all tadpole graphs, and the tree- level 
tadpole coupling A* will never appear explicitly in results; 
its only effect is to precisely cancel the sum of the loop 
tadpoles, which therefore do not need to be calculated 
either. 

There is a subtlety in the case of spontaneous break- 
ing of gauge symmetries, however. The first method 
above applies for arbitrary £, and independence of the 
pole masses with respect to £ provides a valuable check. 
But, if we choose the second method, expanding around 
a VEV that is not the minimum of the tree-level poten- 
tial, then the results obtained below without including 
tadpole graphs are only valid in Landau gauge, £ = 0. 
This is because for £ ^ 0, there is a non-trivial tree-level 
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FIG. 1: One- loop self-energy topologies. 



mixing between the longitudinal components of the mas- 
sive gauge field and the Goldstone bosons (found among 
the Ri), so that the propagator does not have the simple 
form of equation (|2.9[) . This is no great loss in practice, 
since the general two-loop corrections to the effective po- 
tential are only simple and have only been evaluated in 
Landau gauge, anyway. 

In this paper, I will not include tadpole loop graphs, as 
I favor the method of expanding around the VEVs that 
minimize the full effective potential, rather than VEVs 
that minimize only the tree-level potential. This means 
that Landau gauge £ = should be used, whenever it 
makes a difference, in practical applications of this work 
involving massive vector bosons. 3 

The one-loop and two-loop contributions to the self- 
energy function are calculated below as the sum of con- 
nected, one-particle irreducible Feynman graphs with two 
external scalar lines. The topologies for one-loop and 
two-loop self-energy graphs are shown in Figures 1 and 
2. The nomenclature used in this paper for the Feynman 
diagrams and their associated loop functions is as follows. 
Letters A, B are assigned to the one-loop topologies, and 
W, X, Y, Z, S, U, V, M for the two-loop topologies, with- 
out regard to the types of propagators. For each such 
two-loop topology, a canonical ordering is assigned to the 
internal propagator lines, as shown. (No canonical order- 
ing is needed for the one-loop topologies or for the two- 
loop topology S.) Then the loop diagrams are assigned 
a name given by the letter corresponding to the topol- 
ogy, with subscripts S, F, F, or V according to whether 
the internal lines (in the canonical ordering) correspond 
to scalar, chirality-preserving fermion, chirality-violating 
fermion, or vector propagators. Diagrams related by in- 
terchange of the external scalar lines, or by symmetries 
of the topology with respect to interchanges of internal 
lines, are not considered separately. A one-loop example 
and a two-loop example of this naming scheme are shown 
in Figure 3. 

For each Feynman diagram, the integrations over in- 
ternal momenta are regulated by continuing to d = 4 — 2e 



3 However, the calculations of the individual graphs do not actually 
use any particular minimization condition for the VEVs. So one 
could also use the results given here, together with a further 
calculation of the relevant tadpole graphs, to obtain the answers 
valid for an expansion about the VEVs that minimize the tree- 
level potential. 
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FIG. 2: Two-loop self-energy topologies. The numbers corre- 
spond to a canonical ordering for the propagator lines, which 
applies both for subscripts labeling propagator types (as S, 
F, F, or V) and for squared mass arguments in the corre- 
sponding loop integral functions. 





B S v(x,y) 



M V p SF p(x,y,z,u,v) 



FIG. 3: Examples of the nomenclature for self-energy dia- 
gram loop integral functions used in this paper. Dashed lines 
represent scalars, wavy lines vectors, solid lines with a single 
arrow are chirality-preserving fermion propagators, and solid 
lines with clashing arrows are chirality-violating fermion prop- 
agators. The arguments x, y, z, it, v denote the corresponding 
propagator squared masses. 



dimensions, according to 

d 4 k -» (27rp) 2e / d d k. 



(2.14) 



In the dimensional rcgularization scheme, the vector 
bosons also have d components, while in the dimensional 
reduction scheme they have d ordinary components and 
2e additional components known as epsilon scalars, which 
in general have a different squared mass. Counterterms 
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for the one-loop sub-divergences and the remaining two- 
loop divergences are added, according to the rules of min- 
imal subtraction, to give finite loop functions, which then 
depend on the renormalization scale Q given by 



47re~V. 



(2.15) 



The resulting renormalization schemes are known as MS 
|64| and DR respectively. A further coupling con- 
stant redefinition in the latter case removes the depen- 
dence on the unphysical epsilon scalar masses, resulting 
in the DR scheme (§5 appropriate for (softly-broken) 
supersymmetric models. In this paper, the symbol 



1 for MS 
for DR' 



(2.16) 



will be used to present results simultaneously in the MS 
and DR schemes. 

Below, the result for each Feynman diagram is given in 
terms of a corresponding loop integral function, for exam- 
ple Bsv(x, y) and M V p SF p(x, y, z, u, v) for the Feynman 
diagrams shown in Figure 3. [These functions will be 
given explicitly in equations (|3.7|) and (|4.86|l for massive 
vectors, and in equations (|3.10() and l|5.26[) for massless 
vectors.] The variables x, y, z, u, v here stand for squared 
masses. Each loop function associated with a Feynman 
diagram also depends on the external momentum invari- 
ant s, as defined in equation (|1.2|) . and on the renormal- 
ization scale Q, but these are not listed explicitly among 
the arguments, since they are the same for all loop func- 
tions in a given calculation. By convention, the usual 
symmetry factors associated with the Feynman diagram 
are factored out of the corresponding loop function, but 
minus signs associated with fermion loops and other fac- 
tors associated with the evaluation of the diagram are 
not. All counterterms are included within the loop inte- 
gral functions. 

Tarasov has described an algorithm |5l| for reducing 
any two-loop self-energy integral with arbitrary masses to 
a minimal set of basis integral functions. This algorithm 
has been conveniently implemented in the computer pro- 
gram TARCER [5^. In this paper, I use the Tarasov 
algorithm to reduce the scalar two-loop self-energy inte- 
gral functions to a (slightly modified) basis of functions 
that were defined in |62j| . 

Briefly, the basis functions can be described as fol- 
lows. At one loop, integrals with scalar propagators and 
the topologies of A and B in Figure 1 give rise to basis 
functions A(x), B(x,y), which are essentially Passarino- 
Veltman |72( functions. At two loops, integrals with 
scalar propagators and the topologies S, U, and M in Fig- 
ure 2 give rise to basis functions S(x,y,z), U{x,y,z,u) 
and M(x,y, z,u,v). Another basis function T(x,y,z) is 
obtained by differentiating — S(x, y, z) with respect to x. 
Also, I(x, y, z) is defined as the two-loop vacuum integral 
function obtained by setting s = in S(x,y,z). These 
basis functions include one-loop and two-loop counter- 
terms, and they are finite and independent of e. Their 



precise definitions were given in ref. [62|. That paper 
also gave a complete and efficient algorithm for their nu- 
merical calculation, following the strategy put forward in 
^^j-|6]| using the differential equations method of |55|- 

Because of the existence of graphs with the topology 
V shown in Figure 2 with the 2,3 propagators having the 
same squared mass, it is useful as a matter of notational 
convenience to define the function 

d 

V(x,y,z,u) = -—U(x,y,z,u). (2.17) 
dy 

This is not an independent basis function; instead, it can 
be expressed algebraicall y in terms of the basis functions, 
using equation (3.22) of |62l |. 

It is useful to note that all of the functions mentioned 
in the previous two paragraphs have smooth limits as 
squared-mass arguments are taken to zero, with two ex- 
ceptions; there are logarithmic singularities as x — > in 
the basis function T(x, y, z) and as y — ► in V(x, y, z, u). 
These infrared singularities must cancel from physical 
quantities, but as a book-keeping device it is useful to de- 
fine modified versions of these functions that have smooth 
massless limits. Accordingly, in reference 62], the func- 
tion 

T(0, y, z) = lim [T(x, y, z) + B(y, z)Ex] (2.18) 

x— >0 

was defined. Here we have used the notation that loga- 
rithms of dimensionful quantities are always expressed in 
terms of 



hxX = ln(A/Q 2 ). 



(2.19) 



An analytic expression for T(0, y, z) was given in equa- 
tion (6.18) of |62(. Similarly, it is useful to introduce 



V(x, 0, z, u) 



lim 



V{x,y,z,u) 



I(y',z,u) 



, (2.20) 



where the prime indicates a derivative with respect to 
that squared-mass argument. The function V{x 1 Q, z,u) 
can be expressed algebraically in terms of the basis func- 
tions; the result is in Appendix A of the present paper. 
Also, below I will use 



B(x',y) 
I(x',y,z) 



d_ 

dx 
d_ 

dx 



B(x,y) 
I(x,y,z) 



(2.21) 
(2.22) 



to express a few results more compactly than would 
otherwise be possible. These derivatives, expressed al- 
gebraically in terms of the basis functions, are given 
in equations (3.1) and (5.3) of ref. .@2|; they are also 
easy to express analytically in terms of logarithms and 
dilogarithms, using equations (6.3)-(6.9) of that paper. 
Note that B(x,y) and I(x,y,z) are symmetric under 
interchange of arguments, so B(x',y) — B(y,x') and 
I(x',y,z) = I(y,x',z) etc. 
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To summarize, here is the procedure for evaluating the 
two-loop scalar self-energy function for a general renor- 
malizable theory, given the running parameters as in- 
puts. The first step is to minimize the two-loop effec- 
tive potential |63| to find the VEVs of the theory. Then, 
one expands the real scalar, two-component fermion, and 
real vector fields about those VEVs, and diagonalizes the 
tree-level squared mass matrices. This specifies the in- 
teraction Lagrangian in the form of equations H2.lM2.4j l. 
defining the tree-level couplings. These can be plugged 
in to the formulas given in section III, IV, and V below, 
which involve loop diagram functions [like B$v (x, y) and 
Myj: SF -p(x, y, z, u, v)], which are in turn written in terms 
of basis functions: 

A(x), B(x,y), I(x,y,z), S(x,y,z), T(x,y,z), 
T(0,y,z), U{x,y,z,u), M{x,y,z,u,v), (2.23) 



and some additional functions 

B(x',y), I(x',y,z), V(x,y,z,u), V(x,0,z,u) (2.24) 

that can be written algebraically in terms of the basis 
functions. Finally, the basis functions in general are eval- 
uated by computer, using the methods described in detail 
in 



III. ONE-LOOP SCALAR SELF-ENERGY 
FUNCTIONS 

The one-loop self-energy function matrix for scalars in 
the general theory specified in section ITll is: 



ng 5 = ^ kk A s (ml) + lx ikl X^B S s(ml,m^ + Re[y KLi y KLj ]BMrn%,ml) 

+Re[ y K ^y K ' L 'nd KK ,M LLI ]B TT (ml,ml)+g alk g^ k Bsv(mlml) + ±g^A v (m 2 a ) 
+i fl «V M Bw(m2,mg) ) 

where repeated indices are summed over, and the loop integ ral functions arc given by: 



(3.1) 



A s (x) 
Rss(x,y) 
B FF (x,y) 
Bjrp(a;,y) 

A v (x) 



A(x), (3.2) 

-B(x,y), (3.3) 

(x + y-s)B(x,y)-A(x)-A(y), (3.4) 

2B(x,y), (3.5) 



4A(x) + 2xS m -C x [xA{x)], 



(3.6) 



{2x - y + 2s)B{x, y) + A(x) - 2A(y) 



+C y [(x + y- s)A(y) - (x - s) 2 B(x, y)], (3.7) 



B V v(x,y) 



\B{x,y) 



2& 



MS 



~£ x [xB(x,y)] 



+ -C y [yB(x,y)} + -C x C y {xA{y) + yA{x) 

+ [2s(x + y)-x 2 -y 2 - s 2 ]B{x, y)} . (3.8) 

The last function Byv(x,y) combines the contributions 
from the diagram with two vector propagators and the 
corresponding diagram with two ghost propagators. In 
each case, the vector masses are taken to be non-zero. In 
the limit of massless gauge vector bosons, one obtains: 

A v (0) = 0, (3.9) 
B sv (x,0) = (3-S)(x + a)B(0,x) 

+(3-2£)A(x)+2(£-l)a. (3.10) 

As explained in sectional the MS scheme results are ob- 
tained by setting fr^r — 1 in the above equations, while 



the DR scheme results have & 



MS 



0. Note that even 



though the function Bsy(x,y) involves a vector propa- 
gator, it is the same in both schemes, because the vector 
indices are contracted with momenta rather than a met- 
ric tensor. The gauge-dependent parts (acted on by an 
C x or C y ) do not differ between the two schemes, for the 
same reason. 

When computing a pole mass in the two-loop approx- 
imation, one needs ITW^s) = dH^/ds. This is easily 
obtained from the preceding, either by using the analyt- 
ical expression for the function B(x,y), or by using the 
identity: 

a dB fo V) = -xBtf, y) - yB(x, y') - 1, (3.11) 

which follows from equation ljA.2jl and dimensional anal- 
ysis. 



IV. TWO-LOOP SCALAR SELF-ENERGY 
FUNCTIONS 

In this section, I present the results for two-loop con- 
tributions to the scalar self-energy. They are divided into 
subsections, depending on whether the Feynman diagram 
contains internal propagators lines that are only scalars; 
scalars and fermions; scalars and one vector; fermions 
and one vector; or scalars, fermions, and one vector. 
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Throughout the following, repeated indices are summed to the self-energy: 
over. 



A. Diagrams with only scalar propagators 

The two-loop Feynman diagrams that involve only 
scalar propagator lines yield the following contribution 



n g) = ^ kl X kmn X lmn W S sss(mlmlm 2 m ,m 2 n ) + jX^ kl X klmm X SS s(mlmf,m 2 m ) 
+ lx lkl V km X lmnn Yssss(ml mlm 2 m , m 2 n ) + ±X M X mn X klmn Z ssss (mlm 2 ,m 2 m , m 2 n ) 
+ ^\ lklm V klm S S ss(ml,m?,m 2 m ) + ± (A lfei A jfemn + A^A lfcm ") X lmn V S sss(m 2 k , mf, m 2 m , m 2 n ) 
+ lx ikl ^ km X l ^X mn PVsssss(mlmlm 2 m ,m 2 n ,m^ (4.1) 
Here the loop integral functions are given by: 





= [I(x,z,u) - I(y,z,u)]/(y- x), 


(4.2) 




= A(z)[A(x)-A(y)]/(x-y), 


(4.3) 


Yssss(x,y,z,u) 


= A(u)[B(x,z)-B(x,y)]/(y-z), 


(4.4) 


Zssss(x,y,z,u) 


= B(x,y)B(z,u), 


(4.5) 


Ssss(x,y,z) 


= -S(x,y,z), 


(4.6) 


Vssss(x,y,z,u) 


= U(x,y,z,u), 


(4.7) 


Vsssss(x,y,z,u,v) 


= [U(x,y,u,v) -U(x,z,u,v)]/(y- z), 


(4.8) 


M S ssss(x,y,z,u,v) 


= -M(x,y,z,u,v). 


(4.9) 



In several of these functions, denominators threaten to vanish when masses become degenerate. Taking the corre- 
sponding limits, one obtains: 

Wssss(x,x, z,u) = -I(x',z,u), (4-10) 

X sss {x,x,z) = A(z)Ex, (4.11) 

Yssss(x,y,y,u) = -A(u)B(x,y'), (4.12) 

~Vsssss(x,y,y,u,v) = -V(x,y,u,v). (4.13) 



B. Diagrams with scalar and fermion propagators 

In this subsection, I present the results for Feynman diagrams that involve both scalar and fermion propagators, 
but no vector propagators. 

The contributions from diagrams with the topology W in hgurc[21 in which propagators 1,2 are scalars and propa- 
gators 3,4 are fermions are: 

ng } = iA^ fe 'Re[y MWfe / f ' w ' / M M M'M^]W S5 ^(m 2 fe ,m^m2 f ,m2 r ) 

+ l^ kl y MNk yMNiWssFF(m 2 ,m 2 ,m 2 M ,m 2 N ), (4.14) 

where 

W S sff( x > 2/> z , u) = -2W S sss(x,y,z,u), (4.15) 

W S SFF(x,y,z,u) = {(z + u-x)I(x,z,u)-A(x)[A(z)+A(u)]}/(x-y) + (x^y), (4.16) 



8 



Wssff{x, x, z, u) = (z + u — x)I(x' ,z,u) — I(x,z,u) — [A(u) + A(z)]lnx. (4-17) 
The contributions from diagrams of the topology M, taking propagators 1,2,3,4 to be fermions and 5 to be scalar, are: 

ng } = Re[y KMi y LNj y K ' L ' p y M ' N ' p M KK ,M LL ,M MM ,M NN '] M TWF5 (m^, m 2 L ,m 2 M ,m 2 N , m 2 p ) 
+2Rc[y KMl y LNj y KUp y M ' Np M LL ' M M m>] M FFFFS {m 2 K ,m 2 L , m 2 M ,m 2 N , m 2 p ) 
+Rc[(y KMl y LNj + y KMl y LNl ) y K u P yMN> P M LL ' M NN '} M FFFFS {m 2 K , m 2 L ,m 2 M ,m 2 N , m 2 p ) 
+2Re[y KMl y LN iy KLp y M ' N 'PM M M>M NN ,]M FFFFS (m 2 K , m 2 L , m 2 M , m 2 N ,m 2 p ) 

+Rc[y KMl y LN:i yKLpyMNp] M FFFFS {m 2 K ,m 2 L , m 2 M ,m 2 N , m 2 ), (4.18) 

where 

M FFFFS (x,y,z,u,v) = 2M(x,y,z,u,v), (4.19) 

M FFFFS (x,y,z,u,v) = (y + z - v - s)M(x,y,z,u,v) - U(x,z,u,v) -U(u,y,x,v) + B(x,z)B(y,u), (4.20) 

M FFFFS (x,y,z,u,v) = (x + z- s)M(x,y,z,u,v)-U(y,u,z,v)-U(u,y,x,v), (4.21) 

M FFF -Fs(x,y,z,u,v) = {x + y -v)M(x,y,z,u,v) -U(x,z,u,v) - U(y,u,z,v) + B(x,z)B(y,u), (4.22) 
Mffffs(^, y, z, u, v) = (xu + yz — vs)M(x, y, z, u, v) — xU{z, x, y, v) — zU{x, z, u, v) 

—uU(y, u, z, v) — yU(u, y, x, v) + S(x, u, v) + S(y, z, v) + sB(x, z)B(y, u). (4.23) 

The results from diagrams of the topology M, taking propagators 1,3 to be scalars and 2,4,5 to be fermions, are: 

n (2) = X i km ^ [y^j y VPk y N'P'm MLLiMNNiMpp ^ M __ {m 2^ m | ; m ^ ^ mp) 

+2Rc[y LN] y LPk y N ' Pm M NN ,] M SFSFF (m 2 k ,m 2 L ,m 2 m , m 2 N , m 2 P ) 

+Rc[y LNl y L p k y N p, m M pp ']M SFSFF (m 2 k ,m 2 L ,m 2 ml m 2 N ,m 2 P )j + (i <-» j), (4.24) 

where 

M s ^ 5 prp(x,y, z, u,v) = 2M(x,y,z,u,v), (4.25) 

^■Sfsff( x > Vi z ' u ' v ) = ( v ~~ x + v)M{x, y, z, u, v) + U(y, u, z, v) — U(x, z, u, v) — B(x, z)B(y, u), (4.26) 

M SFSFF-( x >y> z > u > v ) = {y + u-s)M(x,y,z,u,v)-U(x,z,u,v)-U(z,x,y,v). (4.27) 

The contributions from diagrams of topology V, with lines 1,2,3 taken to be scalar and 4,5 to be fermions, are: 
n (2) = x iU \i km {Re[y NPl y N ' p ' m M NN ,M Ppi ]W SSSFF {m 2 ^ 

+Re[y NPl y NPm \V S ssFF(m 2 k ,m 2 , m 2 m , m 2 N ,m 2 P f) , (4.28) 

where 

V SSSFF (x,y,z,u,v) = -2V S ssss(x,y,z,u,v), (4.29) 

VsssFF(x,y,z,u,v) = {(y-u-v)U(x,y,u,v) + [A(u)+A(v)]B(x,y)}/(y-z) + (y<-*z), (4.30) 

VsssFF(x,y,y,u,v) = {u + v - y)V(x,y,u,v) + U(x,y,u,v) + [A(u) + A(v)]B(x,y'). (4.31) 

The results from diagrams of topology V, with lines 1,2,3,4 taken to be fermions, and 5 to be scalar, are: 

ng> = 2Rc[y KLi y K ' M ^y L ' N Py M ' N 'PM KK ,M LL ,M MM/ M NN ,]V Fms ( m 2 K , m 2 L , m 2 M ,m 2 N , m 2 p ) 

+2Re[(y KLl y K ' M i + y KL >y K ' Mt )y L N P y M ' Np M KK ,M MM ,]V FFFFS (m 2 K , m 2 Ll m 2 Ml m 2 N , m 2 p ) 
+2Re [y KL l y K ' M > y L Np y M N> P M KK , M NN ']V FFFFS (m 2 K , m| , m 2 M , m 2 N , m 2 p ) 
+2Re[y KLl y KM] y L ' Np yM>N P M LLI M MM '] V FFFFS (m 2 K , m|, m 2 M , m 2 N , m 2 p ) 
+2Rc[{y KLl y KM3 + y KL3 y KMl ) yLN P yM>N> P M MM ' M NN ']V FFFFS {m 2 K ,m 2 L ,m 2 M7 m 2 N ,m 2 p ) 
+2Rc[y KLl y KMl y LNp y MN P]V FFFFS (m 2 K ,m 2 L ,m 2 M , m 2 N ,m 2 p ), (4.32) 
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where 

V FFFF S (x,y,z,u,v) = -2V S ssss(x,y,z,u,v) 7 (4.33) 
Vffffs( x >V> z > u > v ) = {( v -y- u )U(x,y,u,v) + [A(v)~A(u)]B(x,y)}/(y-z) + (y^z), (4.34) 
^FFFFs( x 'y^ z ^ u ' v ) = 2[zU(x,z,u,v) ~ yll(x,y,u,v)]/(y - z), (4.35) 
V F FFFs( x 'y' z ' u > v ) = {{y + u - v )[{ s - x -y)U{x,y,u,v) + I{y,u,v)] 

+ [A(u) - A(v)][(s -x- y)B(x,y) + A(y)}} /2y(y - z) + (y «-» z) 
+{(u- v)[2xT(x,u,v) + 2S(x,u,v) - A{x) - A{u) - A(v) +x + u + v- s/4] 
+u(x + u — v — s)T(u,v,x) — v(x + v — u — s)T(v,x, u)}/2yz, (4.36) 
V ' FFTFs( x ^y^ z ^ u ^ v ) = [( s - x - y)U(x,y,u,v) + I(y,u,v)]/(y - z) + (y ^ z), (4.37) 
Vffffs{x, V, z, u,v) = {(y + u- v) [(s - x - y)U(x, y, u, v) + I(y, u, v)] + yS (x, u, v) 

+ [A(u) - A(v)][(s -x- y)B{x,y) + A(y)]}/2(y - z) + (y <-> z). (4.38) 

In the case y = z, the necessary limits of these functions are given by: 

y ffffs( x > y> V> u > v ) = (y + u - v)V(x, y, u, v) - U(x, y, u, v) + [A(v) - A(u)]B(x, y'), (4.39) 

y ffffs( x > y> y> u > v ) = 2 y y ( x i v> u > v ) - 2U ( x , y, u > v), ( 4 - 4 °) 

Y FFFFs( x ^y'y^ u ' v ) = {y(y + u ~ v)[(x + y- s)V(x,y,u,v) +I(y',u,v)] + [(u-v)(x - s) - y 2 ]U(x,y,u,v) 
+ (u - v)[2xT(x,u,v) + 2S(x,u,v) - I(y,u,v) - A(x) - A(u) - A(v) + x + u + v - s/4] 
+u(x + u — v — s)T(u, v, x) — v(x + v — u — s)T(v, x, u) 

+[A(u) - A(v)][(s -x- y)yB(x,y') + (x - s)B(x,y) + y}}/2y 2 , (4.41) 
Y FFTFs( x ^V^y^ u ^ v ) = {x + y- s)V(x,y,u,v) - U(x,y,u,v) + I(y' ,u,v) (4.42) 
Vffffs(x, y,y,u,v) = {{y + u- v) [(x + y- s)V(x, y, u, v) + I(y', it, v)] + (s - x - 2y - u + v)U(x, y, u, v) 

+S(x, u, v) + I(y, u, v) + [A{u) - A(v)] [(s - x - y)B(x, y') - B(x, y) + Ey]}/2. (4.43) 

The limit y — * appropriate for massless fermions is only needed when the corresponding propagator has no mass 
insertions. For the case V FFFFS (x, y, y, u, v), this limit is trivial, since yV(x, y, u, v) vanishes as y — > 0. The remaining 
case is given by: 

Vffffs(%, 0, 0, it, v) = { (it — v)(x — s)V(x, 0, u, v) + (s — x — u + v)U (x, 0, u, v) + S(x, u,v) + 1(0, u, v)} /2 

+ [A(u) - A(v)][sB(0, x') + A{x)/x\. (4.44) 

This function is well-defined as s — > x, although (x — s)V(x,0,u,v) and B(0,x') are both singular in that limit. 



C. Diagrams with scalar and one vector propagators 

The contributions from two-loop Feynman diagrams with scalar propagators, exactly one vector propagator, and 
no fermion propagators are 

Ilg> - ^ kl 9 akm g alm Wsssv(mlmf,m 2 m ,ml) + \^ kl g aakl X S sv (ml , mf, m 2 a ) 

+ \x kl V km 9 aalm ^sssv(ml mlm 2 m , m 2 a ) + l -g alk g a ' l \ klmm Y vsss (ml m 2 k , m 2 ,m 2 m ) 

+ 1 {g aik X jkmn + g ajk yikmn ^rnn^^^^^^j 

+\ ikl \ jkm 9 aln 9 amn Vssssv(mlml m 2 m ,m 2 n , m 2 a ) + ^ k g^ l X kmn X lmn Vvssss(m 2 a ,mlml m 2 m ,m 2 n ) 
+ 1 {g aik x jkl + g°*\M)g°™\i™V S vsss(ml,ml, m 2 , m 2 m , m 2 n ) 

+ ^ x ikrn xj ln g akl g a m n Mssssv{m 2 jm ^ m 2 mj ^^2) 

+ {g ail X jkm + g a 3 l ykra )g akn X l ™M VS SSS (™l , m\ , mf , , m£ ) , (4.45) 
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where the loop functions are given by 

W S ssv{x,y,z,u) = {(2x + 2z-u)I(x,z,u) + A(x)[A(z)-2A(u)] + 2xA(u) 

+C U [-(x - z) 2 I(x, z, u) + A(u){(z + u- x)A(x) + xA(z)}] }/(x - y) + (x <-+ y), (4.46) 

Wsssy (x, x, z, u) = (2x + 2z- u)I(x', z, u) + 2I(x, z, u) + [A(z) - 2A(u)] Ex + 2A(u) 

+£„{ — (x — z) 2 I{x ', z, u) + 2(z — x)I(x, z, u) + A(u)[(z + u — 2x)lnx + x + A(z)]}, (4.47) 

X SS v(x,y,z) = (4A(z) + 2z5 m -C z [zA(z)])[A(x)-A(y)]/(x-y), (4.48) 

X SS v(x,x,z) = {4A(z) + 2zS ms -jC z [zA(z)])]^x, (4.49) 

Ysssv (x, y, z, u) = {4A(u) + 2u6 MS - C u \uA{u)\) [B(x, z) - B(x, y)]/{y - z), (4.50) 

Y S ssv(x,y,y,u) = - {4A{u) + 2u<5 M g- C u [uA{u)\) B(x,y'), (4.51) 
Yvsss(x,y,z,u) = A{u){(2y - x + 2s)B(x,y) + A(y) + C x [yA{x) - (y - s) 2 B(x,y)]} /(y - z) + (y <-» z), (4.52) 
Y VSS s{x,y,y,u) = A(u){(2y - x + 2s)B(x, y') + 2B(x, y) + hry 

+C X [A(x) -(y- s) 2 B(x, y') - 2(y - a)B(x, y)] }, (4.53) 
Usvss{x,y, z,u) = {(z - u)[(x - s)U(x,y, z,u) + 2xT(x 1 z 1 u) + 2S(x,z,u) - I(y,z,u) - A(x) - A(z) 
— A(u) + x + z + u — s/4] + z(x + z — u — s)T(z, u, x) — u(x + u — z — s)T(u, x, z) 
+ [A(z) - A(u)][(x - s)B(x, y) - A(y)]}/y + C y {(z - u)[(s - x)U(x, y, z, u) + I(y, z, u)} 

+ [A(z) - A(u)][(s - x)B{x, y) + A(y)]}, (4.54) 
Vssssv(x,y,z,u,v) = |(d - 2y- 2u)U(x,y,u,v) + [2A(v) - A(u)]B(x,y) 

+C V [(y - u) 2 U(x 1 y, u,v) - yS(x,u, v) + (y - u - v)A(v)B(x, y)] |/(y - z) + (y <-» z), (4.55) 

Vssssv{x,y,y,u,v) = (2y + 2u- v)V(x,y,u,v) - 2U{x,y,u,v) + [2A(v) - A(u)]B(x,y') 

+C v {2(y -u)U {x, y, u, v) - (y - u) 2 V(x, y, u, v) - S{x, u, v) + A(v) [{y-u- v)B(x, y') + B{x, y)] } , (4.56) 

Vvssss{x,y,z,u,v) = {(x - 2y - 2s)U(x,y,u,v) - I(y,u,v) 

+C X [(y - s) 2 U {x, y, u, v) - yS{x, u, v)] }/(y - z) + (y <-► z), (4.57) 

Vvssss{x,y,y,u,v) = (2y - x + 2s)V(x,y,u,v) - 2U(x,y,u,v) - I(y' ,u,v) 

+C X [2{y - s)U{x, y, u, v) - (y - s) 2 V{x, y, u, v) - S{x, u, v)] , (4.58) 

Vsvsss(x,y,z,u,v) = |(u - v)[(x - s)U(x,y,u,v) - I(y,u,v)\ + [A(u) - A(v)][(x - s)B{x,y) - A(y)]j/y(y - z) 
+ (y <-> z) + j(v - u)[2xT(x,u,v) + 2S(x,u,v) - A(x) - A(u) - A(v) + x + u + v - s/4] 
+u(s + v — x — u)T(u, x, v) + v(x + v — u — s)T(v, x, u)^/yz + C y |^|(u — v)[(s — x)U(x, y, u, v) 
+I(y,u,v)] + [A(u) - A(v)][(s - x)B(x,y) + A(y)]}/(y - z) + (y <-> z)] , (4.59) 

V|y L sss (x, y, y, u,v) = |(t> — u)[{x — s)U {x, y, u, v) — I(y, u, v) + 2xT(x, u, v) + 2S(x, u, v) — A{x) — A{u) 
— A(v) + x + u + v — s/4] + u(s + v — x — u)T(u, x, v) + v(x + v — u — s)T(v, x, u) 
+ [A(u) - A(v)][{s - x)B{x,y) + A{y)]}/y 2 + J 1 £„{( M - v)[(x - s)U(x,y,u,v) - I{y,u,v)] 

+ [A(u) - A(v)][(x - s)B(x, y) - A(y)]}, (4.60) 

Y svsss( x > V' V> u ' v ) = {( v ~ u )[y{x - s)V(x, y, u, v) + (x - s)U(x, y, u, v) + yl(y', u, v) - I(y, u, v) 
+2xT(x, u, v) + 2S(x, u, v) - A(x) - A(u) - A(v) + x + u + v - s/4] + u(s + v - x - u)T(u, x, v) 
+v(x + v-u- s)T(v,x,u) + [A(u) - A(v)][(x - s)yB(x,y') - (x - s)B{x,y) -y]}/y 2 , (4.61) 

M VS sss(x,y, z, u, v) = (s-x + y + z-2u + v)M(x,y, z,u,v) -U(x, z,u,v) + 2U(z,x,y,v) 
-U(u,y,x,v) + U{y,u,z,v) - B(x, z)B(y,u) + C x (s - z)[(v - y)M(x,y, z,u,v) 
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+U (x, z, u, v) - B(x, z)B(y, u)] + (v - y)U(u, y, x, v) + S(x, u, v) - A(x)B(y, u) , 
M SS ssv{x, y, z, u, v) = ((x - v/4 - s/2)M(x, y, z, u,v)-U (x, z, u, v) + B(x, z)B{y, u)/A 



(4.62) 



(y - x)[zM(x, y, z, u, v)/2 + U(z, x, y, v)\ + S(x, u, v)/2 - A(v)B(x, z)/2 



+(x <-> y and z <-> u) + (x <-> z and y <-> u) + (x <-> u and y <-> z). 



(4.63) 



D. Diagrams with fermion and one vector propagators 

This subsection contains the contributions to the scalar self-energy coming from two-loop diagrams that involve 
fermion propagators, one vector propagator, and zero scalar propagators. 

The results for diagrams with topology M with propagators 1,2,3,4 taken to be fcrmions and 5 to be vector are: 



ng } = Re[y KM % LNj g a «' g a N r M KK ,M LL ' M MM ,M 



aK' „oM ' 



rLL' - 



rNN'~ 



M- 



FFFFV 



(m 2 K ,m 2 L ,m 2 M ,m 2 N ,m 2 a ) 



+2Re[y KMi y LN ig a K L ' g a N M ' M LL/ M MM ,]M FTWFV (m 2 K , m\, m 2 M ,m 2 N ,m 2 a ) 
+Re[{y KMi y LN i + y KMj y Lm )g a K L ' <fi?' ' M LL ,M NN ,] M FFFFV (m 2 K ,m 2 L , m 2 Ml m 2 Nl m 2 a ) 
+2Re[y KMl y LN3 g a K L g%¥'MMM'M NN ']M FFFFV (m 2 K ,m 2 L ,m 2 M , m 2 N ,m 2 a ) 
+Re[y KMt y LNj gfg a ^]M FFFFV {m 2 K , m 2 L ,m 2 M , m 2 N ,m 2 a ), 



where 



FFFFV 



(x,y,z,u,v) = 8M(x,y, z,u,v) - £ v [2vM(x,y, z,u,v)], 



(4.64) 
(4.65) 



M- 

M FWFy (x,y,z,u,w) = |(2y - v - s)M(x,y, z,u,v) - 2U(x, z, u, v) + B(x, z)B(y, u) +5^ 

+C V [{v(s/2 - x) + (x - y)(z - u)/2}M(x, y, z, u, v) + (z - u)U(x, z, u, v) + (v + x - y)U(z, x, y, v) 
-S(x,u,v)/2- S(y,z,v)/2 + A(v)B(x,z)]} + (x «uandy <-> z), (4.66) 

M FFFFV (x,y,z,u,v) = {(4x-2s)M(x,y,z,u,v) - 4U(y,u,z,v) + [2B(y,u) + l]<%jg 

+C V [v(s/2 — x)M(x, y, z, u, v) + vll(y, u, z, v)] j + (x <-> z and y <-» u), (4.67) 

m ffTfv( x > V> z > u > v ) = ( 2x ~ v)M(x, y, z, u, v) - 2U (x, z, u, v) + B(x, z)B(y, u) + 5 m 

+C v ^2zx(x — y — v)M(x, y, z, u, v) + [s(u — z — v) + (x + z)(v — u) + z(3x — 2y + z)]U(x, z, u, v) 

+2x(v — u)T(x, u, v) + u(s + v — x — u)T(u, x, v) + v(x + v — u — s)T(v, x, u) 

+ (2v -2u- z)S(x,u, v) + (u — z — v)I(u, z, v) + A(v)[(x + z - s)B(x, z) - A(z) +u — v] 

+v[v + x - A{x) - A(u) - s/A]^/2z + (x <-> y and z ^ u), (4.68) 

M FFFFV (x, y, z,u, v) = ^(s - x - z)[(s - y - u)M(x, y, z,u,v)/2 + 2U(x, z, u,v)] + S(x, u,v) + I(x,y, v) 
+[(x + z - s)B(x, z) - 2A{x) + 2x + v/2 - s/4]6 m 

+C V [sz(x — y) + x(2xu — 2yz — vu)]M(x, y, z, u, v) + [s(z — u — v) + (x + z){u + v) 
+z(x — 2y — z)]U(x, z, u, v) + {x(s — x — 3?i + 3v)T(x, u, v) + v(v — s)T(v, x, u) 

+ (v - 4x + 3y - s)S(x,u,v) + A(v)[x - A(x) - v] + v[v + x - 2A(x) - 3I(x,y,v) - s/4]}/3 /2^j 
+(x <-» y and z <-> u) + (x <-> z and y <-> u) + (x <-> u and y <-> z). (4.69) 
The results for diagrams of topology V with lines 1,2,3,4 taken to be fermion and 5 to be vector are: 



n 



(2) 



2Re 
+2Re 



yKLiyK'Mj g aL' g aM' MkkiMll/Mmm/M 



NN' 



v 7 



„KLi„,K' Mi „aN „aN' 

y y " " 



( 2. 2 2 2 2\ 

FFFFV\ m K> m Li m Mi m Ni m a) 

9T 9 a M m kk' M NN , V FF FFV (m% , m\ , m 2 M ,m 2 N ,m 2 a ) 
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+2Rc [(y KLi y KMj + V KL] Vkm 1 ) g a L N g a ^' M mm ' M nn ^ V FFFFV (m 2 K , m 2 L ,m 2 M , m 2 N , m 2 a ) 
+2Rc ^yKLiyK'Mj +y KL jy K'Mi^ g aN g aM> Mkk ,M mm ,} V FFFFV (m 2 K , m|, m 2 M , m 2 N , ml) 
+2Re [y KLl y K Mjg a N L ' g a A ^M LL ,M MM '] Y FFFFV (m 2 K , m\, m 2 M , m 2 N , m 2 a ) 

+2Rc [y KLl y K M 3 g a L N g a N M ] V F FFFv(m 2 K ,m 2 L ,m 2 M ,m 2 N ,m 2 a ), (4.70) 

where the integral functions are (including cases with degenerate fermion masses y — z): 

VFT?FTv( x >y> z > u > v ) = {-&U(x,y,u,v) + 4B(x,y)& m + £ v [2vU(x,y,u,v)]} /(y - z) + (y <-» z), (4.71) 

y ffffv ( x > 2A 2A U ' v ) = 8V ( X > 2/, u ' w ) + 4B ( X > 2/') (5 ms ~ C v [2vV(x, y, u, v)}, (4.72) 

^FFFFv( x ^y^ z ^ u ^ v ) = 2y{-4U(x,y,u,v) + [2B(x,y) + 1]6 MS + £ v [vU(x,y,u,v)]}/(y - z) + (y <-» z), (4.73) 
^FFFFvi x ^y^y^ u ^ v ) = 8yV(x,y,u,v) - 8U(x,y,u,v) + [4yB(x,y') + 4B(x,y) + 2}8 m 

-£ v {2v[yV{x, y, u, v) - U(x, y, u, «)]}, (4.74) 
VFFFTv( x >y' z ' u ' v ) = l 4 ( s ~ x ~ v)U{x,y,u,v) +4I(y,u,v) + 2[(x + y- s)B(x,y) - A(y) + y]S ME 

+£ v [v(x + y - s)U (x, y, u, v) - vl(y, u, v)]}/(y - z) + (y <-> z), (4.75) 
^f f~ffv i x ' 2/' U ' w ) = 4 K X + 2/ - s)V(a;, y, u,v)-U (x, y, u, v) + I(y, u, v)} + 2[(x + y- s)B(x, y') 

+B(x,y) - A(y)/y]6 m + £ v {v[(s -x- y)V(x,y,u,v) + U(x,y,u,v) - I(y',u,v)]} , (4.76) 



V ffffv ( x > 2/' z ' u ' v )= 2 ( v -y- u)U(x, y, u, v) + 2[A(v) - A(u)]B(x, y) + y[2B(x, y) + 1]S- 



MS 



+£ v {[v(y + u)-{y- u) 2 ]U{x, y, u, v) + yS{x, u, v) + (u - y)A{v)B{x, y)}] /{y -z) + {y^z), (4.77) 
v ;ffffv( x >2/'2A u ^) = 2 (y + u- v)V(x,y,u,v) - 2U(x,y,u,v) + 2[A(v) - A(u)]B(x,y') 

+[2yB(x, y 1 ) + 2B(x, y) + l]S m + £ v { [(y - u) 2 - v(y + u)]V{x, y, u, v) + (2u - 2y + v)U(x, y, u, v) 

+S(x, u, v) + (u - y)A(v)B(x, y 1 ) - A(v)B(x, y) } , (4.78) 
VFFFFv( x 'y> z > u ' v ) =({y + u- v )[( s ~ x ~ y)U(x,y,u,v) + I(y,u,v)] 

+ [A(u) - A(v) - y&Mg] [(s - x - y)B(x, y) + A(y)]j /y(y - z) + {y <-► z) + (v(s + u - x - v)T{y, x, u) 

+u(x + u-v- s)T(u,x,v) + (u- v)[2xT(x,u,v) + 2S(x,u,v) - A(x) - A(u) - A(v) +x + u + v- s/4]j/yz 

+ s ms + c v{(l v (v + u ) - (y~ u f]\( x + y - s)U(x,y,u,v) - I(u,v,y)} + {y- u)A(v)[(s -x- y)B(x,y) 

+A(y)]j /2y(y — z) + (y <-» z) + (uv(s + u — v — x)T(v, x, u) + u 2 (x + u — s — v)T(u, x, v) 

+u(u - v)[2xT(x,u,v) - A(v)} + (yz + 2u 2 - 2uv)S(x,u,v) + uv[A(x) + A(u) -x-v + s/4]^)/2yz|, (4.79) 

V ffffv( x > 2A2A«^) = (y(y + u - v)[(x + y - s)V(x,y,u,v) + I(y',u,v)] + [(u-v)(x - s) - y 2 ]U(x,y,u,v) 
+(v - u)I(y, u, v) + \A(u) - A(v)] [y(s - x - y)B(x, y') + (x - s)B(x, y) +y}+ v(s + u - x - v)T(v, x, u) 
+u(x + u-v- s)T(u,x,v) + (u- v)[2xT(x,u,v) + 2S(x,u,v) - A(x) - A(u) - A(v) +x + u + v- s/4]) jy 2 
+[{x + y- s)B(x, y') + B{x, y) - A(y)/y]S m 

+£ v {y[v(y + u) -(y- u) 2 ][(s -x- y)V{x,y,u,v) - I(y',u,v)} + [s{y 2 + uv - u 2 ) 

+y 2 {2u + v - x - 2y) + ux(u - v)]U(x,y, u,v) + (y 2 + uv - u 2 )I(y,u,v) + A(v){y(y - u)(s - x - y)B(x, y') 
+ [u(s — x) — y 2 ]B(x, y) + y[A(y) + y — u]} + uv(s + u — v — x)T{v, x,u) + u 2 (x + u — s — v)T(u, x, v) 

+u(u- v)[2xT(x,u,v) - A(v)] + (y 2 + 2u 2 - 2uv)S(x, u, v) + uv[A(x) + A{u) -x-v + s/4]|/2y 2 , (4.80) 
V ffffv {x, y, z, u, v)= {[y + u- v) [(s - x - y)U (x, y, u, v) + I(y, u, v)] 

+[A(u) - A(v) - y5 m ] [(s-x- y)B{x, y) + A{y)fj /(y - z) + (y <-» z) + S(x, u, v) - 2A{v) 
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+[x- A(x) + y + z + 2u- 3s/4]5 m + £ v {(jv(y + u) - (y - uf][(x + y - s)U(x,y,u,v) - I(y,u,v)} 

+ (y - u)A(v)[(s -x- y)B{x, y) + A{y)^J /2(y - z) + (y z) + (v(3x -3u-v + s)T(v, x, u) 

+u(3x -u-3v + s)T(u, x, v) + 2x(x - s)T{x, u, v) + (8x + 3y + 3z - 7u - 4v - 2s)S(x, u, v) 

+ [2A{x) - AA{u) -2x + u + v]A(v) + v[A(x) + A{u) +x-2u-v + s/4]j /6|, (4.81) 

^ffffv{x, y, y, u, v ) = (y + u- v)[(x + y- s)V(x, y, u, v) + I(y', u, v)) + (v - x - 2y - u + s)U(x, y, u, v) 
+I{y, u, v) + \A(u) - A(v)] [(s - x - y)B(x, y') - B(x, y) + hiy] + S(x, u, v) - 2A(v) 
+[y(x + y- s)B(x, y') + (x + 2y- s)B(x, y) - A(x) - 2A(y) +x + y + 2u- 3^/4]^ 

+C v ^[v(y + u) - (y - u) 2 ][(s - x - y)V(x, y,u, v) - I(y', u,v)] + [s(2y - 2u - v) + uv + 2ux + vx 

+Auy + 2vy - 2xy - u 2 - 3y 2 ]U(x, y,u, v) + [2y - 2u - v]I(y,u, v) + A(v)[(y - u)(s - x - y)B(x 1 y') 

+(s — x — 2y + u)B(x, y) — y + (2y — u)lny] + (v(3x — 3u — v + s)T(v, x, u) + u(3x — u — 3v + s)T(u, x, v) 

+2x(x - s)T(x, u, v) + (8x + 6y- 7u - 4v - 2s)S(x, u, v) + [2A(x) - 4A(u) - 2x + u + v]A(v) 

+v[A(x) + A(u) + x - 2u - v + s/4]) /sj /2. (4.82) 

In the last case, the limit y = z = needed for massless fermions is slightly non-trivial: 

Vffffv( x i 0; 0, u, v) = (u — v)(x — s)V(x, 0, it, v) + (s — x — u + v)U(x, 0, u, v) + S(x, u, v) + 1(0, u, v) 
+2[A(u) - A(v)][sB(0,x') + A(x)/x] - 2A(v) + [{x - s)B(0,x) - A(x) + x + 2u - 3s/4]^jg 

+C v ^u(u — v)(x — s)V(x, 0, u, v) + [uv + 2ux + vx — u 2 — 2us — vs]U(x, 0, u, v) — (2u + v)I(0, u, v) 

+A(v) [(s + 2u- x)B(0, x) - 2uxB(0, x')] + (v(3x - 3u - v + s)T(v, x, u) + u(3x - u - 3v + s)T(u, x, v) 

+2x(x - s)T(x, u, v) + (8x -7u-4v- 2s)S(x, u, v) + [2A{x) - AA(u) - 2x - 5u + v]A(v) 

+v[A(x) + A(u) + x-2u-v + s/4]j /sj /2. (4.83) 

This function is well-defined as s — > x, even though B(0, x') and (x — s)V(x, 0, u, v) are each singular in that limit. 

E. Diagrams with scalar, fermion, and vector propagators 

In this subsection I present the results for two-loop diagrams that contain one vector propagator and both fermion 
and scalar propagators. 

The contributions of diagrams of topology M with line 1 taken to be a vector propagator, lines 2,4,5 to be fermion, 
and line 3 to be scalar, are: 

ng } = 2{lm[g aa y KM > g a N K ' y M ' m M KK ,M M M>]M v7s7F {m 2 a ,m 2 K ,mlm 2 M ,m 2 N ) 
+Im[g ail y KM 3g a N K 'y M N>iM KK ,M NN '] U vTsFT (m 2 al m 2 K , m 2 , m 2 M , m 2 N ) 
+Im[g aa y KM >g a K N y M ' N ' l M M M'M NN/ ] M VFsTF (m 2 a , m 2 K , mf, m 2 M ,m 2 N ) 
+Im[g aa y KM ^g]fy Mm ]M VFSFF (m 2 a ,m 2 K ,m 2 ,m 2 M ,m 2 N )) + (i ~ j), (4.84) 



where 



M vfsff( x ' V> z,u,v) = (y - x + 2z - 2u + v)M(x, y, z, u, v) - B(x, z)B(y, u) - U(x, z, u, v) 

+2U (z, x, y,v) + U (y, u, z, v) - 2U(u, y, x,v)+£ x (x - y + v)[(s - z)M(x, y, z, u,v) + U (u, y, x, v)} 

+(s-z)[U(x,z,u,v)-B(x,z)B(y,u)] + S(x,u,v)-A(x)B(y,u) , (4.85) 
m vfsff( x ' y,z,u,v) = (y- z-2u + v- s)M(x, y, z, u, v) - U(x, z, u, v) + 2U (z, x, y,v) + U (u, y, x, v) 
-B(x,z)B(y,u) +£ x (z - s)[(y + z - v - s)M(x,y, z,u,v) -U(x, z,u,v) + B(x,z)B(y,u)] 
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+{[s(x + 3y — v) + uv — ux — uy + vy — xy — y 2 — 2yz]U(u, y, x, v) + x(s — x — u + v)T(x, u, v) 
+v(u + v — x — s)T(v, x, u) + 2u(v — x)T(u, v, x) + (y — 2x + 2v)S(x, u, v) + (x + y — v)I(x, y, v) 

+ (s -u- y)A(x)B(y, u) + A(x)[A(y) + x - v] + x[A{u) + A(v) -x-u + s/4]}/2y , (4.86) 

Mypg-pp^x, y, z, u, v ) = (x — y + 2u — v — 2s)M(x, y, z, u,v) + U (x, z, u, v) — 2U(z 1 x, y,v) — U (y, u, z, v) 

+B(x, z)B(y, u) + C x (x + y - v)[(s - z)M(x, y, z,u, v) + U(u, y, x,v)] + (z - s)[U(x, z, u,v) 

-B{x, z)B{y, «)] - S(x, u, v) + A(x)B(y, u)] , (4.87) 

MypsFF(xj y, z, u i v ) — [s(2z — 2u — v) + (2u — v + 2x — y — 2z)u — yz]M(x, y, z, u, v) + (u + z)U (x, z, u, v) 
+(y + 2u — 2s)U (u, y,x,v) — 2uU (y, u, z, v) + (it + s)B(x, z)B{y, u) — S(x, u, v) + {[s(x + y — v) 
—2ux — x 2 + v z + xz — yz]U(z, x, y, v) + 2(v — y)zT(z, y, v) + y(s + v — y — z)T(y 1 z, v) 
+v(v + z-s- y)T{v, y, z) + (x - 2y + 2v)S(y, z, v) + [A(v) - A(y)][(z - s)B(x, z) - A{x)\ 
+(y - x - v)I(x,y,v) + (y - v)[A{y) + A(z) + A(v) - y- z-v + s/A]}/x 

+C X (z - s)[(yz -yu + vu- vs)M(x, y, z,u, v) + (u - z)U(x, z, u,v) + (y - v)U(z,x,y, v) 

+(s - u)B(x, z)B(y, u)] + [s(v + x + y) + u(y - x - v) + y(y - x - 2z - v)]U(u,y,x, v)/2 

+ [x(x + 3v — 3u — s)T(x, u, v) + v(3x + v — 3u — s)T(v, x, u) + 2u(s — u)T(u, x, v) 

+ {Ax - 3y + 6z - 8u + Av - As)S(x,u,v)]/6 + (y + u - s)A{x)B{y,u)/2 

+(z - s)[A(y) - A(v)]B(x, z) + A(x)[2A(v)/3 - A(u)/3 - A(y)/2] + (x - y + v)I{x,y, v)/2 

+ [(2u - x - v)A{x) + x{2v + x - u - s/4 - A(u) - A(v)}]/6 . (4.88) 

The contributions from diagrams of topology V with line 1 taken to be vector, lines 2,3 to be scalar, and lines 4,5 
to be fermion are: 

+Rc[y MNk y M Ni] VvssFF(ml, m 2 kl ml m 2 M , m 2 N )j , (4.89) 

where 

y vss¥f( x > V' z ' u ' v ) = ~ 2 Vvssss{x, y, z, u, v), (4.90) 
Vvssff(x, y,z,u,v) = | (y - u- v)[(x -2y- 2s)U(x, y, u, v) - I(y, u, v)] + 2yS{x, u, v) 

+[A(u) + A(v)][(x - 2y - 2s)B(x,y) - A(y)]}/(y - z) + (y <- z) 

+C 

X {(y-s) 2 

{{y — u — v)U(x, y,u, v) + [A(u) + A(v)]B(x,y)} / (y - z) + (y <-» z) 

+{u(u — 3x + 3v — s)T(u, x, v) + v(3u — 3x + v — s)T(v, x, u) + 2x(s — x)T(x, u, v) + (5s — 5a; — 3y 

— 3z + 4u + 4v)S(x, u, v ) + 2[A{u) + A(v)][x - A{x)\ + {2x-u- v)A(x) + x(s/2 -2x-u- u)}/3 j, (4.91) 

VvssFF{x,y,y,u,v) = (y - u - v)[(2y -x + 2s)V(x,y,u,v) - I(y',u,v)] + [x- Ay + 2u + 2v- 2s]U(x,y,u,v) 
-I{y, u, v) + 2S{x, u, v) + [A{u) + A{v)} [{x -2y- 2s)B{x, y') - 2B{x, y) - Ey] 

+C X { (y - s) 2 {(u + v- y)V{x, y, u, v) + [A{u) + A(v)]B(x, y')} + (y - s)(3y - 2u -2v- s)U (x, y, u, v) 

+2(y - s)[A{u) + A(v)]B(x, y) + {u(u - 3x + 3v - s)T(u, x, v) + v(3u - 3x + v - s)T(v, x, u) 
+2x(s - x)T(x, u, v) + (5s - 5x - 3y - 3z + 4u + Av)S(x, u, v) + 2[A(u) + A(v)} [x ~ A{x)} 

+ (2x — u — v)A{x) + x(s/2 -2x-u- u)}/3|. (4.92) 

Finally, the results for diagrams of topology V with lines 1,3 taken to be scalars, 2 to be vector, and 4,5 to be 
fermion, are: 

njf = 2Im[( 3 « fe A^+.^ fe A* fc ^^ (4-93) 
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where 

V SVSFF -{x,y,z,u,v) = ((y + u-v)[(s - x)U(x,y,u,v) + I(y,u,v)} + [A(u) - A(v)][(s - x)B(x,y) 
+A(y)]j /y(y — z) + (y «-> z) + (u(x + u — v — s)T(u, x, v) + v(s + u — x — v)T(v, x, u) 
+ (u - v)[2xT(x,u,v) + 2S(x,u,v) - A(x) - A(u) - A(v) + x + u + v - s/A]j/yz 
+£ y { ({v-y- u)[(s - x)U(x, y, u, v) + I(y, u, v)} + [A(v) - A(u)][{s - x)B(x, y) + A(y)fj /{y - z) 

+fo <->*)}. (4-94) 
^y SF p{x,y,y,u,v) = ((v-y-u)[(s - x)U{x,y,u,v) +I(y,u,v)\ + [A(v) - A(u)][(s - x)B{x,y) + A(y)} 
+u(x + u — v — s)T(u, x, v) + v(s + u — x — v)T(v, x, u) + (u — v)[2xT{x, u, v) + 2S{x, u, v) 
-A(x) - A(u) - A(v) +x + u + v ~ s/4]j /y 2 

+ -^^C y [(y + u - v)[{s - x)U(x, y, u, v) + I(y, u, v)] + \A(u) - A(v)][(s - x)B(x, y) + A(y)]}, (4.95) 

V svWf( x ' y ' y ' u ' ^ = { y ( y + u ~ v )[( x ~ s ) v ( x > Vi u , v ) + Hv', u i v)]+u(x + u-v- s)T(u, x, v) 

+v(s + u — x — v)T(v, x, u) + (u — v)[(x — s)U {x, y, u, v) — I(y, it, v) + 2xT(x, u, v) + 2S(x, u, v) — A(x) 
-A(u) - A(v) +x + u + v- s/4] + [A(u) - A(v)][(s - x)yB(x,y') + (x - s)B(x,y) + y]\/y 2 . (4.96) 



V. TWO-LOOP FUNCTIONS INVOLVING MASSLESS VECTORS 

In the preceding formulas, the vector bosons have been taken to have generic non-zero masses. However, when 
the corresponding gauge symmetry is unbroken, the vector boson will be massless. Also, the three-point couplings 
of such massless vector bosons occur only when the other two particles are degenerate in mass. In this section I will 
give the appropriate limits of the two-loop integral functions for massless vector bosons. It should be noted that the 
presentations of some of the results in this section are not unique, because of the existence of identities involving basis 
functions with vanishing squared masses [see equations HA.llfl - ljA.20fl in the Appendix]. 

For the integral functions of subsection IIV CI one finds: 



W S ssv(x,x,x,0) = 37(0, x, x) - A{x) + 2x, (5.1) 

X S sv(x,x,0) = 0, (5.2) 

Ysssv(x,y,y,0) = 0, (5.3) 

Yvsss(.0,y,z,u) = A(u){(3 - £)(y + s)B(0, y) + (3 - 2£)A(y)}/(y - z) + (y <-> z), (5.4) 

Y VSS s(0,y,y,u) = A(u){(3-0[(y + s)B(0,y') + B(0,y)] + (3-2^y}, (5.5) 

Vsvss{x,0, z,z) = 0, (5.6) 

Vssssv(x, y, y, y, 0) = AyV(x, y, y, 0) - 2U(x, y, y, 0) - A{y)B{x, y') + (1 - £)T(0, x, y), (5.7) 

Vvssss(0 : y,z,u 7 v) = [(£ - 3)(y + s)U(0, y,u,v) + (2£ - 3)I(y,u, v)]/{y - z) + (y «-> z), (5.8) 

Vvssss(0,y,y,u 7 v) = (3-OKy + s)V(0,y,u,v)-U(0,y,u,v)} + (2^-3)I(y',u,v)}, (5.9) 

V S ysss(^,0, z,u, u) = 0, (5.10) 
M S sssv(x, x, y, y, 0) = 2(x + y - s)M(x, x, y, y, 0) - 2U(x, y, y, 0) - 2U(y, x, x, 0) + B(x, y) 2 

+(2£-2)T(0,x,y), (5.11) 
Myssss(0, y, z, u, y) = (2y + z-2u + s)M(0, y, z, u,y) + U (y, u, y,z)-U (u, y, 0, y) - U(0, z, y, u) 

+2U(z, 0, y, y) - B(0, z)B(y, u ) + (l-£){(z + s) [U(0, z, y, u) - B(0, z)B(y, u)} 

-2yT(y, 0, u) - 2uT(u, 0, y) - 45(0, y, u) + 21{y, z, u) + 2[s - A{z)]B{y, u) 

+2A(u) + 2A(y) - 2u - 2y + s/2}/(z - s). (5.12) 



The massless- vector limits of the integrals of subsection IIV Dl arc: 



Mp FFFv (x,x,y,y,0) = (6 + 2£)M(x,x,y, y, 0), 



(5.13) 
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M 



FFFFV 



[x,x,y,y,0) 



M 



FFFFV 



M 



FFFFV 



{x,x,y,y,0) 



(x,x,y,y,0) 



M F p F p V (x,x,y,y,0) 



(1 + Ol(x +y- s)M(x, x, y, y,0)-U [x, y, y,0)-U (y, x, x, 0)] 
+2(1 - OT(0, x, y) + 2B(x, yf + 2S m , 

(3 + £)[(x + y- S )M {x, x, y, y,0)-U {x, y, y,0)-U {y, x, x, 0)] 
+ [4B(x,y) + 2]5 ws , 

2x(l + OM(x, x, y, y, 0) + [(1 - £)(x + y - s)/y - 4}U(x, y, y, 0) + 2B(x, y) 2 
+ (1 - i)[AxT{x, 0, y) + 2yT{y, 0,x) + (x + y- s)T(0, x, y) + 45(0, x, y) 
-1(0, y, y) - 2A(x) - AA{y) + 3(x + y - s/2)]/y + 2S m , 
2[(x + y-s) 2 + (^- l)xy]M(x, x, y, y, 0) + (3 + - x - y)[U(x, y, y, 0) 
+U(y, x,x,0)]+2(l-£)[xT(x,0,y)+ yT(y, 0,x)- A(x) - A(y) +x + y- a/2] 
+ (6 - 2£)5(0, x, y) + (l+ 0[I(O, x, x) + 1(0, y, y)} 
+4[(x + y- s)B(x, y) - A(x) - A(y) +x + y- s/A]S m , 



and 



V 



FFFFV 



(x,y,y,y,0) = (6 + 2£)V{x, y, y, 0) + 4B(x, y')S m , 



Vffffv( x > V> V> y> 0) = (6 + 2^ [yV(x, y, y, 0) - U(x, y, y, 0)] + [4yB(x, y') + 4B{x, y) + 2}S m , 
y FFFFV (x, y, y, y, 0) - (3 + £) [{x + y- s)V(x, y, y, 0) - U(x, y, y, 0) - A{yf /2y 2 } 
+2[{x + y- s)B{x, y') + B[x, y) - A{y)/y]S m , 



V- 



FFFFV 



(x, y, y, y, 0) - 2£T(0, x, y) + 2[£A(y) - (1 + £)y]B(x, y') - (1 + + 2B(x, y)] 



+ {2yB(x,y') + 2B(x,y) + l]S m , 
Vp-pppvix, y, y, y, 0) = £{ (x + y - s)T(0, x, y) + (s-x + y)T(y, 0, x) + 2xT{x, 0, y) + 25(0, x, y) - A(x) 
+2x + 2y- 5s/4} + (x + y - a)[£A(y) - (1 + y') + [£(s - x)A(y)/ 2/ - (1 + £)y]B{x, y) 



(5.14) 
(5.15) 

(5.16) 

(5.17) 

(5.18) 
(5.19) 

(5.20) 

(5.21) 



+ (1 - U)A{y)\ /y+[(x + y- s)B(x, y') + B(x, y) - A(y)/y}S m , 
Vffffv{x, y, y, y,0) = £[{x + y- s)T(0, x, y) + 2yT(y, 0, x) + 5(0, x, y) + I{0, y, y)] 



(5.22) 



+(x + y - 8)[£A(y) - (1 + Z)y\B(x, y') + [£A(y) + (1 + £)(«- a: - 2y)]B(x, y) + (1 - £)A{x) 

+ (2 - £)A(y) -x-3y+(3 + £)s/4 + [y{x + y- s)B{x, y') + (x + 2y - s)B{x, y) - A{x) 

-2A{y) + x + 3y- 3s/4]S m . (5.23) 

Finally, the integrals introduced in subsection IIV El have the following forms for massless vectors: 
m vfsff(°> V> z ' u > V) = 2 (V + z ~ u)M{0, y, z, u,y) + U (y, u, y,z)-U (0, z, y, u) + 2U (z, 0, y, y) 
-2U(u, y, 0, y) - B(0, z)B(y, u) + (1 - [(s - z)M (0, y, z, u, y) + U(u, y, 0, y) 
+ {(z + s)U(0, z, y, u) - 2yT(y, 0, u) - 2uT(u, 0, y) - 45(0, y, u) + 21{y, z, u) 

-(z + s)B{0, z)B(y, u) + 2[s - A(z)]B(y, u) + 2A(y) + 2A(u) - 2u - 2y + s/2} /{z - s)l , (5.24) 

M vfsTf(°' z > u ' V) = - M vfsff(°> V' z > u > y) + 2 £K Z _ s)M{0, y, z, u,y)-U (u, y, 0, y)}, (5.25) 
M vfsff(°' y> z > u > V) = ( 2 V - z - 2u ~ S ) M (0, y, z, u,y)-U (0, z, y, u) + 2U {z, 0, y, y) + U(u, y, 0, y) 
-B(0, z)B(y, u) + (1 - [S(y, y, z) - uU(y, u, y, z) + [A(u) - A{y) - z/2 - s/2]B(0, z) 

-A(z)/2 + \{uz + ys - zs)U(0, z, y, u) + {s-y- u)[yT(y, 0, u) + uT(u, 0, y) + 25(0, y, u) 

+{ A(z) - s}B(y, u) - A{u)} + {y-z + u)I(y, z, u) + {sz -su- yz)B(0, z)B{y, u) 

+ {y + u- s)A{y) - 5s 2 /8 + (10u + lOy + 3z)s/8 - (y + uf}/{z - s) /y, (5.26) 

M VFS ff(0, y, z, u, y) = [s(2z - y - 2u) - yz + 2u{u -y- z)]M(0, y, z, u, y) - 2uU (y, u, y, z) 

+ (u + z)U{0, z, y, u) + (s -2u + z)U (z, 0, y, y) + (y + 2u - 2s)U(u, y, 0, y) + S{y, y, z) - 5(0, y, u) 

-1(0, y, y) + (u + s)B(0, z)B(y, u) + (1 - £) [s(y, y, z) - uU(y, u, y, z) - A(z)/2 

+ [A(u) - A(y) - z/2 - s/2]B{0, z) + {(z 2 - yz - su)U(0, z,y,u) + (y + u- s)[yT{y, 0, u) 
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+uT(u, 0, y) + {A(z) - s}B(y, u) - A(u)] + (2u + 2y - z - s)S(0, y, u) + (z - y - u)I(y, z, u) 
+ (zu + ys- s 2 )B(0, z)B(y, u) + (s - y - u)A{y) + 3s 2 /8 - s(10y + lOu + z)/8 



+ (y + u) 2 }/(z- s) 



(5.27) 



and 



Vvssff(0, y, z, u, v) = [(£ - 3) (y + s) [{y-u- v)U(0, y, u, v) + {A(u) + A(v)}B(0, y)] 
+ (2£ - 3)[(y - u - v)I(y, u, v) + {A(u) + A(v)}A(y)]\ /(y - z) + {y <- z) 

+(3-£)S(0,u,v) + (l-Os/2, (5.28) 
VvssFF{0 7 y,y 7 u,v) = (3 - £){(y + s)(y - u - v)V(0, y, u, v) + (u + v-2y- s)U(0,y,u,v) + S(0,u,v) 
- [A(u) + A(v)} [(y + s)B(0, y') + B(0, y)] } + (2£ - 3) [(y-u- v)I(y', u, v) + I{y, u, v) 
+{A(u) + A(v)}hiy] + (1 - Os/2, (5.29) 

Vsvsf~f( x > z, u, u) — £ (s — x)[U(x, z,u,u) — U(x,0,u,u)] + I(z,u,u) — 1(0, u,u) jz. (5.30) 

Certain combinations of the above functions involving massless vector bosons often appear together, and will 
be noted for future reference. In the self-energy functions of scalars that are neutral under the gauge symmetry 
corresponding to the massless vector boson, the combinations 

Gss(x, y) = Mssssv (x, x, y, y, 0) + Vssssv (x, y, y, y, 0) + Vssssv(y, x, x, x, 0) 

= 2(x + y- s)M(x, x, y, y, 0) + 4yV(x, y, y, 0) + 4xV(y, x, x, 0) - 4U (x, y, y, 0) 

-4U(y, x, x, 0) - A(y)B(x, y') - A(x)B(x',y) + B(x, yf , (5.31) 

VFFFFv(x,y,y,y,0) + yV F pp FV (x,y,y,y,Q) - 2yV FFFFV (x, y, y, y, 0) 
-M FFFFV (x, x, y, y, 0)/2 + yM FFFTv (x, x, y, y, 0) - xyM FFFF V (x, x, y, y, 0)/2 + (x «-> y) 

-(x + y- s) 2 M(x, x, y, y, 0) + 4(x + y - s) [U (x, y, y, 0) - yV(x, y, y, 0)] + S(0, x, y) 

+ (3x + y- s)T(x, 0, y) + (x + y)B(x, y) 2 + (y - x - s)A(x)B(x, y)/x 
+2(s -x- y)A(x)B(y, x') - 5/(0, x, x) + 6A(x) - Ax - s/4 

+S m {2y(s -x- y)B(x, y') + (s - 4x)B(x, y) + 4A(x) + 2x - s/4}] + (x y), (5.32) 
2V FFFFV (x 1 y, y, y, 0) - V FFFFV (x, y, y, y, 0) - yV TWFFv (x, y, y, y, 0) 
+M FFFFV (x,x,y,y,0) - M FFFFV (x,x,y,y,0) + (x <-> y) 
= (2s - 4x)M(x, x, y, y, 0) + 8U(x, y, y, 0) - 8yV(x, y, y, 0) + 2B(x, yf - 4A(x)B(y, x') 

+S m {-4B(x,y)-4yB(x,y')}]+(x^y) (5.33) 

will appear. They are each independent of the gauge parameter £. For scalars that transform non-trivially under the 
gauge symmetry corresponding to the massless vector boson, the following combinations also occur: 



G F F(x,y) 



G FF (x, y) 



Gssss(x,y,z,u) 
G S sFF(x,y,z,u) 



G 



SSFF 



(x,y,z,u) 



Vssssv(z, u, u, u, 0) + Vvssss(0, x, y, z, u) + My S sss(0, u, x, z, u) 
+M vssss (0,u,y,z,u) 7 

Vy S SFF(0, x, y, z, u) + V FFFFV (z, u, u, u, 0) + uV FFFFV (z, u, u, u, 0) 
-2uV FF - FFV (z,u,u,u, 0) + Mvfsff(0,u, x, z, u) + M V fsff(0,u, y, z,u) 
-!iM„-f C i7-f(0, u, x, z, u) — uM, 
V,, 



l vfsff(°' m ' x ' z ' u ) - uM vfsff(°' M ' V> Z ' U )' 



VSSFF (0,x,y,z,u) - V FFFFV (z,u,u,u,0) + 2V FFFFV (z,u,u,u,0) 
u, u, u, 0) + M yFS -p7p(0, u, x, z, u) + M yFS -pp(0, u, y, z, u 



FFFFV^ ' ' ' ' / ' VFSFF 

- m vfsff(°' u ' x - z ' u ) ~ m vfsff(°' u > V' z > u )- 
These are not gauge-invariant, but satisfy: 

d d d 

-^G ssss (x,x,y,z) = -B S s(y,2) — — B S y(x,0) + 



(5.34) 

(5.35) 

(5.36) 
(5.37) 
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^G S sFF{x,x,y, z) 



0_ 



d d 

-B FF (y, z) — —B sv (x, 0) 
d d 

-B TF (y,z) — —B S v(x,0) 



(5.38) 
(5.39) 



where, from equation (|3.1U|) and the explicit forms of A(x) and B(0,x) given in equations (2.10) and (6.4) of [62 



— — B sv {x, 0) = 2 + 2E(x - a) - Ear • 



(5.40) 



The ellipses in each case refer to terms that vanish as s — > x. These identities are useful for checking gauge invariance 
of the pole mass. ^ 



VI. A SUPERSYMMETRIC EXAMPLE 

In this section, I will work through the details of an 
example intended to serve both as a consistency check 
and as a point of reference for the preceding results. 

Consider a supersymmetric model with three chiral su- 
perfields $o> and with charges 0, +1,— 1 under 
a U(l) gauge symmetry. The complex scalar compo- 
nents of the supcrficlds can be written in terms of real 
scalars i?o,+,- and Zo,+,- as 00 = (-Ro + iiio)/v2 and 
<t>+ = (R+ + H+)/V2 and 0_ = (R- + UJ)/y/2, and 
their two-component Weyl fermionic components are ipo, 
ip + , and tp_ respectively. Also, there is a U(l) gauge 
field A^ , and a gaugino two-component Weyl fermion \- 



These fields couple to the charged fields with strength g. 
The superpotcntial is given by: 



(6.1) 



In 

M±. Then 



Let us take Mo, M±, and Y to be real and positive 
the following, I will denote x = Mq and z 
the real scalars Ro, Iq, R+, I+, J_ have squared 
masses x, x, z, z, z, z; the Weyl fermions ipo, 4>+, ip-, X 
have masses y/x, y/z, y/z, and 0; and the gauge field A 
is massless. The interaction Lagrangian density is given 
by: 



-C = rVs((M-05 + c - c -) + *V*(fo + 0o)(l^+| 2 + I'M 2 ) + Y 2 (l<M 2 |<M 2 + |0+I 2 |0o| 2 + |<M 2 |0o| 2 ) 

W{\<t>+? - \4>- 1 2 ) 2 /2 + Y [(0oMM + MM^o + <t>-Mo) + c.c] + V2g [(</>; X M ~ <f>-X>l>-) + C -C-] 
+igA»{<t> + dr<l>% - r+d^ + + <f>*_d^ - j-d^t) + 5 2 (|<M 2 + |0-| 2 )A"^ 

- 5 ^[^f^+-^vM- (6.2) 
This determines the couplings of the theory as specified in eqs. (I2.1I) - (I2.4[I . for example: 

X R + R + R-R. =Y 2_ X R R + R- = Y y/^/2, y++Xl + = _ ig> gAR+I + = _ g AI + R + = g> ( g 3) 

etc. 

Applying the results of section ITTTI to this model, one obtains for the one-loop self-energy functions of the neutral 
scalars: 

n£> o = F 2 [2A s (z) + xB SS (z,z) + B FF (z,z) - z%(z,z)], (6.4) 

n 2fl = Y 2 [2k s {z) + {x + Az)B ss {z,z) + B FF {z,z) + zB- ¥¥ (z,z)], (6.5) 

<k = (6.6) 

It follows immediately from equations l|3.3|) and l|3.5|l that the difference between IIjj and B^ qRq vanishes identically. 
Since there is also no mixing in the one-loop self-energy matrix between Ro and Iq, one can write a single self-energy 
function for the neutral complex scalar field. Putting equation (|6.4|l [or (|6.5|l ] in terms of the one-loop basis functions 
gives simply: 

n$ = -Y 2 (x + s)B(z,z). (6.7) 
At two-loop order, we have, by direct application of the results of sections llVI and Ivl 

n / ib = 2Y i S SSs(x,z,z)+Xsss(z,z,x)+Xsss(z,z,z) + (x + 2z)Wssss(z,z,x,z) 
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and 



+2x\Jssss(z, z, x, z) + xY S sss(z, z, z, x) + xY S sss(z, z, z, z) + xZ S sss(z, z, z, z)/2 
+ (x 2 + 2xz)V S ssss(z, z, z, x, z) + xzM S ssss(z, z, z, z, x) + YJ S sff(z, z, x, z) 
+xV S ssff(z, z, z, x, z) - 2zV FFFFS (z, z, z, x, z) - 2zVp F p FS (z, z, z, z, x) + zY FF - FFS (z, z, z, x, z) 
+zY FFFFS (z, z, z, z, x) + Y FFFFS (z, z, z, x, z) + Y FFFFS (z, z, z, z, x) + zM FFFFS (z, z, z, z, x) 
zM FFFF g(z, z, z, z, x) xzYLg-pcjp-p (z, z, z, z, x) -t- xMg F g FF (z, z, z, z, x)] 
+g 2 Y 2 [2X S ss(z, z, z) + 2xY ssss(z, z, z, z) - xZ S sss(z, z, z, z) + AW S sff(z, z, 0, z) 

+AxY S ssff(z, z, z, 0, z) - 8zV FFFFS (z, z, z, 0, z) + 4zV FFFFiS (z, z, z, 0, z) + 4Y FF ffs(z, z, z, 0, z) 
+2W S ssv (z, z, z, 0) + xG S s(z, z) + G FF (z, z) - zGp-p(z, z)], (6.8) 

(21 (21 a 

n floflo " n Wo =^z[AVssss(z,z,x,z) + 2Yssss(z,z,z,x) + 2Yssss(z, z, z, z) + Z S sss(z, z, z, z) 
+ (6x + 4z)V SS sss(2 : , z, z, x, z) + (3a; + 2z)M S ssss(z, z, z, z, x) + 2Y S ssff(z, Z, z, x, z) 
-\-2xY ggg FF (z, z, z, x, z) + 2 V FFFF g (z , z, z, z) + 2V FFFF g (z, z, z, z, a;) 
+M FFFF g,(z, z, z, z, x) 4- xYlgp-gp-p^z , z, z, z, x) + 2Mg F £, FF (z, z, z, z, a;)] 
+4g 2 y 2 z[Z S 555(z, z, z, z) + 2Y S 5 SS (z, z, z, z) + 4V FFFFS (z, z, z, 0, z) + 4Vsssff(-z, z, z, 0, z) 
+4M SFSFF (z, z, z, z, 0) + GfesC*. z) + Gpp(z, z)/2]. (6.9) 

( 2) _____ 

and n^ oFo = 0. The right-hand side of equation (|6.9|l vanishes identically when written in terms of the basis functions. 
Therefore, we again obtain a single expression for the self-energy function of the neutral complex scalar at two-loop 
order. Written in terms of the basis functions, equation (|6.8|l takes the simple form (with = 0): 

II? = Y 4 {xB(z, z) 2 — 2(xz + sz + sx)M(z, z, z, z, x) + (x + s)[2U(z, z, z, a;) — 4zV(z, z, z, a;)]} 

+2.g 2 r 2 (a; + s){(2z - s)M(z, z, z, z, 0) + 2T(0, z, z) + 4{A(z) - z]B(z', z) - 4B(z, z) - 2}. (6.10) 

Let us now consider the self-energy function of the charged scalar fields. At one loop order, the diagonal elements 
are all the same: 

n Sfl + = 4+7+ = n fl!fl_ = n^j_ = Y 2 [A s (x) + A s (z) + (x + 2z)B ss (x 1 z) + B FF (x, z)] 

+ g 2 [A S (z) + 2B FF (0,z) + B S y(z,0)] (6.11) 

while the off-diagonal elements vanish, thanks to supersymmetry: 

^R+R- = -^i- = Y 2 ^x~z[2B SS {x, z) + Bpp(x, z)] = 0. (6.12) 

It follows that the charged scalars have a common self-energy function, 

= ~Y 2 (s + z)B(x, z) + g 2 [4zB(0, z) + (1 - £){(« + z)B(Q, z) + 2A(z) - 2s}}. (6.13) 

(The term proportional to 1 — £ vanishes in the limit s — > z.) The two-loop corrections are, from sections IV and V: 

( 2*1 A 

U r+r + = Y [Ssss(x,x,z) + Ssss(z,z,z) +X S ss(z,z,x) + X S ss(z,z,z) + 2X SS s(x,x, z) 

+ (x + 2z){Wssss(x, x, z, z) + Wssss(z, z, x, z) + Yssss(x, z, z, x) + Y S sss(x, z, z, z) 
+2Yssss(2, ar, a:, z) + Z S sss(x, z, x, z) -I- 2Vssss(x 1 z, x, z) + 2Vssss(z, x, z, z) 
+Vsssff(x, z, z, x, z) + Y S ssff(z, x, x, z, z)} + (a; 2 + 4a;z + 8z 2 )Y S ssss(z, x, x, z, z) 
+ (x 2 + 8xz + -Iz^Ysssssix, z, z, x, z) + (8zx + 4z 2 )M S ss SS (x, z, z, x, z) + W S sff(x, x, z, z) 
+W SSFF (z, z, x, z) -I- zV FFFF5 (x, z, z, x, z) + zV FFFFS (x, z, z, z, x) + 2xV FFFFS (z, x, x, z, z) 
+Y FFFFS (x, z, z, x, z) + Y FFFFS (x, z, z, z, x) + 2Y FFFFS (z, x, x, z, z) + 
-\-2tZ V ' sss~f~f{z •> %t %t %) H~ 2 ; -M-p 1 -prp^i^(x, z, z, x, z) + xM^^r-p^^^z, x, x, z, z) 
~ (~2x]VI^^i^"^r^i ^x 7 z ? z ? z^ "I - 2z]VI^(^i^^t^i (z , x ? z ? z^ — I - 2zlVIg^i gj^pT^x 3 z 7 z 7 x 7 z^j 
+.g 2 Y" 2 [X S s5(^,^,x) + 2X S 5 S (z, z, z) - 2S 5 5 5 (z,z,z) + (4z - 2x)\J S sss{z, x, z, z) 

+ (x + 2z){W S sss(z, z, x, z) + Y S sss(x, z, z, z) + 2V S ssff(x, z, z, 0, z) + G S 55s(z, z, x, z)} 
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+2W S sff(z : z, 0, z) + W S sff(z : z, x, z) + 2zY F - FFFS (0, z, z, x, z) + 2zV F j^ FS (0, z, z, z, x) 
+2zY F - FFFS (x 1 z, z, 0, z) + 2Y FFFFS (0 1 z, z, x, z) + 2Y FFFFS {0, z, z, z, x) + 2Y FFFFS (x 1 z, z, 0, z) 
-4M FFFFS (0, z, z, x, z) + AzM SFSFF (z, 0, x, z, z) + 4zM SFSFF (x, z, z, 0, z) + Y VSS s(0, z, z, x) 
+Yy 555(0, z, z, z) + W S ssv(z, z, z, 0) + G S sff(z, z, x, z)] 
+g 4 [X S ss(z, z, z) + 3Ssss(z, z, z) + 2~W S sff(z, z, 0, z) + AY FFFFS (0, z, z, 0, z) + 8Y FFFFS (z, 0, 0, z, z) 
+4zM F - FFFS (0, z, z, 0, z) + 4zV fTTf5 (0, z, z, 0, z) + W SSS y (z, z, z, 0) + Yy S ss(0, z, z, z) 
+2G SS ff(2,z,0,z) + ...], (6.14) 

with n^ iJ _ and nf 3 7 and IT^ the same. (The ellipses in equation l|6.14[) indicate that the contribution of order 

g is not complete, since in this paper I have not included the contributions with two vector lines. Terms of order <? 4 
are therefore consistently dropped from here on.) Also, 

n i?ji^ = - n /+/_ = 2Y 4 y/xz[W S sss(z,z 1 x 1 z) + 2\J S sss{x,z,x 1 z) + 2\Jssss(z,x 1 z 1 z) +Y S sss(x,z,z,x) 
+Ys555(£, z, z, z) + 2Y S sss(z, x, x, z) + Z SS ss(x, z, x, z) + (x + 4:z){Y S ssss(x, z, z, x, z) 
+Vsssss(z, x, x, z, z) + M S ssss(x, z, z, x, z)} + W SSFF (z, z, x, z)/2 + Y FFFFS (x, z, z, x, z) 
~\~V F ~ F ~ FF S {x, z, z, z, x) -\- 2V FFFF g(z, x, x, z, z 

) + Vsssff(x, z, z, x, z) + Y S ssff{z, x, x, z, z) 
+zV sssippix, z, z, a;, z) 4- zVsssF~F~i z i x > x > z > z ) ^fff~fs( x > z > z > x > z ) ~^ ^-sfsf~f( x > z > z > x > z ) 
~\~^-sfs~ff( z ' x > x > z > z ) "t" z ^s~fsf~f( x > Z > Z > X > z )] 
+g 2 Y 2 ^/xz[2Y S sss(x, z, z, z) - 2Wssss(^, 2, x, z) - W SSF -p(z, z, sc, z) + 4Vsssff(2;, 2, z, 0, z) 
-4M F -p FF - s (0, z, z, x, z) + 4M SFSFF (x, z, z, 0, z) + 4V F - FFF5 (x, z, z, 0, z) 

+2G SS ss(^, z, x, z) + G SSFF (z, z, x, z)]. (6.15) 

Writing these expressions in terms of the two-loop basis functions with &^ — 0, we discover that equation l|6.15() 
vanishes identically, and equation l|6.14(l becomes: 



n (2) _ y4 



|z_B(ir, z) 2 + (s + z)[U(x, z, x, z) — 2zV(x, z,x, z) + U(z,x, z, z) — 2xV(z, x, z, z)\ 

— (zx + z 2 + + 3sz)M(x, z,z,x,z)\ + g 2 Y 2 ^ (8z 2 — 62; z — 2sir + 8sz)M(0, z, z,x, z) + 8z 2 V"(0, z, x, z) 

-4zU(0, z, x, z) + (6z + 2s)[U(z, 0, z, z) - 5(0, z)B(x, z)] + 2(s + z)[T(0, x, z) - 2B(x, z) 

+2{A(z) - z}B(a;, z') - 1] + (1 - £)a} + 0(g 4 ) (6.16) 



where 



A = 2z(s — 3z)V(0, z, x, z) + (s + z)T(0, z) + 45(0, x, z) — 4zl(x, z, z ) — 2I(x, z, z ) -8z 2 5(0,z')5( x, z) 
+ (2s + 6z)B(0, z)B(x, z) + i[A(z) - 2z - s]B(x, z) - 2A(x) - 2A(z) + 2x + 2z - s/2 
+ {(4z 2 + xz + 4zs - sx)U(0, z, x, z) + 2(x 2 - 8z 2 )T(x, 0, z) + 2z(x - 8z)T(z, 0, x) 

+8z[(x - 2z)A{x) - xA{z) - xz]B(0, z)/x - 8z 2 }/(a; - 4z) (6.17) 

embodies the difference between the Landau gauge and Feynman gauge results. (It should be remarked that A can be 
rewritten in various ways that look quite different, using the expressions for 5(0, z') and V(0, z, x, z) in terms of the 
basis functions; see equations (3.1) and (3.22) of |63].) The function A can actually be evaluated entirely analytically, 
using the expressions given in section VI of The limiting expression for small s — z is particularly simple and 
useful: 

A = -2zB(x, z)[2 + 2hi(z - s) - hiz] + . . . . (6.18) 

Here, the ellipses represent terms that vanish as s — > z, which can be consistently neglected since they make only a 
three-loop order contribution to the self-energy function. 

The pole squared masses x p and z p of the neutral and charged scalar fields can now be obtained by 

.,. 1 1 = 1 TT0V * ' 



n ^ 55 x + le^W + (16^ [ u lo (*) + n^(*)n^(*)j , (6.19) 



Z + ^ s X + 16^ n 2 {Z) + (i6^ {Z) + H ^ {Z) \ ■ (6 ' 20) 
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[Compare to equations (|1.4|) - (|1.6|) .] Applying the results above, one finds 
nJHW = -2F 2 £(z,z)x, 



(6.21) 



11^ (x) + IL^(x)Il^'(x) = Y 4 x[~(Az + 2x)M(z, z, z, z, x) + 4U(z, z, z, x) - 8zV(z, z, z, x) + {3B(z, z) 

-8zB(z', z) - 4}B(z, z)] + g 2 Y 2 x[(8z ~ Ax)M (z, z, z, z, 0) + 8T(0, z, z) 
+16{A(z) - z}B(z', z) ~ l6B(z, z) - 8] (6.22) 



with all basis functions on the right evaluated at s = x, and 
I&Hz) = -2Y 2 B(x 1 z)z + 4g 2 B(0,z)z, 



(6.23) 



n^ 2 ] (z) + (z)n^'(z) = Y 4 z [-(42 + 2x)M (x, z, z, x, z) + 2U(x, z, x, z) - 4zV(x, z, x, z) + 2U(z, x, z, z) 

-AxV{z, x, z, z) + {3B(x, z) - AxB{x' ', z) - AzB(x, z') - 4}B(x, z)] + 8g 2 Y 2 z [(2z - x)M(0, z, z, x, z) 
+z{V{0, z, x, z) + B{0, z')B(x, z)} + U(z, 0, z, z) - U(0, z, x, z)/2 + T(0, x, z)/2 

+{A(z) - z + zB(0, z)}B{x, z') + {1 + xB(x', z) - W{x, z)/2}B(0, z) - 1/2] + 0{g 4 ) (6.24) 

with all basis functions evaluated 4 at s — z. These results satisfy several non-trivial consistency checks, as follows. 

First, the mere fact that we can write diagonal self-energy functions H<p and H<p ± is actually a supersymmetric 
consequence of the fact that the real and imaginary components of each of <po, </>+ and 0_ reside in supermultiplets. 
(In a non-supersymmetric theory, the real and imaginary components of the uncharged field 4>o receive different 
self-energy corrections, and there is mixing in the self-energy functions between </> + and 4>*_.) 

Next, consider the fact that the pole squared masses x p and z p must be gauge invariant, since they are physical 
observables. This requires 



T (i)' 



9 ff 



0: 



5 TT 



0. 



(6.25) 



These equations are indeed seen to be satisfied, since £ 
does not appear in equations l|6.21|l - 16.24(l . Let us exam- 
ine how this happened. For the neutral scalars, the one- 
loop gauge invariance is trivial, while the two-loop gauge 
invariance can also be seen from the fact that equation 
(|6.8|l is written in terms of the gauge-invariant combina- 
tions Gss( z , z ), Gff{z,z) and G-p-p(z,z) that were de- 
fined in equations (|5.31() . 15.32(1 and (|5.33() . To see how 
the gauge invariance of the charged scalar pole squared 
mass z p comes about, note that using the analytical ex- 
pressions for A(z) and B(0, z) in equation l|6.13|l gives 



9 n£> (s) = ,g 2 (l - -)[(z + s)ln(z - a) - zlnz], (6.26) 



from which the limit 



9 £-<P± 
follows immediately, and 



a n«(z) = o 



(6.27) 



d_ r 



4>± 



4 As s — * z, the functions V(0, z,x,z) and B(0, z')B(x, z) each 
have a logarithmic divergence, but their sum is well-defined and 
can be evaluated analytically using equations in section VI of 

Ha. 



-2g 2 Y 2 zB(x, z) [2 + 21n(z - s) - lnz] + 



(6.28) 



where the ellipses indicate terms that vanish as s — > z 
and terms of order g A . Combining equation H6.28|l with 
(|6.16l) and l|6.18|l gives the desired smooth limit as s — > z: 



0_ 



n 



(i) 



o. 



(6.29) 



This cancellation explains the absence of terms propor- 
tional to (1 — £) in equation (|6.24|) . 

As another check, suppose that x = 0; then the mass 
of the neutral fermion ip is zero at tree-level, and is pro- 
tected from corrections by a chiral symmetry. It follows 
from supersymmetry that the scalar squared mass also 
vanishes, so x p = 0. This is checked, since equations 
(|6.21|) and l|6.22|l each vanish when x = 0. 

Similarly, suppose that z — 0. Then the masses of the 
charged fermions vanish, and we must have z p = to all 
orders in perturbation theory. Again, this checks, since 
equations (|6.23|1 and l|6.24|l each vanish when z = 0. 

Another important consistency check on the preceding 
results is provided by renormalization group invariance of 
the pole masses of the scalars. The beta functions of the 
parameters of the theory can be written in the general 
form: 



0x = Q 



dX 



1 



dQ 16tt : 



-/3 (1) 
■Px 



1 



(16tt 2 ) : 



-3 {2) 



where 



/?« = 4Y 2 x, 



, (6.30) 



(6.31) 
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[l&g 2 Y 2 - 8Y 4 ]x, 
[AY 2 - 8.g 2 ]z, 
[32g 4 - 8Y 4 ]z, 

3y 3 - 4 5 2 r, 

-6Y 5 + 4g 2 Y 3 + 16g 4 F, 
8.g 5 - 4^ 3 y 2 . 



(6.32) 
(6.33) 
(6.34) 
(6.35) 
(6.36) 
(6.37) 
(6.38) 



The requirement that the pole squared masses must not 
depend on the renormalization scale Q implies: 



d 

dQ' 
d 



(6.39) 
(6.40) 



Let us check these. First, for the neutral scalars, we have 
from equation H6.7(l : 



„2Y 2 x Q ° B [ Z ' Z) - Al3 Y YxB(z, z) 



-2Y 2 (3 X 



dQ 

B(z, z) + x 



dB(z,z) 
ds 



-4Y 2 x{3 z B(z,z'), 



(6.41) 



where all functions on the right side are to be evaluated 
at s = x. Using equations (fP)l . (|6~3TJ| . (|6~33|) . (|6~35|) 
and l|3.11|l . this becomes 

= "4F 2 , + I i ?5 {yM8-20B(z > z)] 
+ 16.g 2 y 2 [B(z,z) + 2z J B(z',z)]}a; (6.42) 



up to terms that contribute only at three-loop order. 
Meanwhile, applying equations 1A.1I) - (|A.10|I . we obtain 



n^WnW'^)] =8r 4 [s(z,z) 

-2zB(z',z) - l]x, 



(6.43) 



Q^n£ } (x) = Y 4 [8 + 12B(z, z) + 16z5(z', z)]x 

-16.g 2 F 2 [l + B(z, z) + 2z5(z', z)]x, (6.44) 

again up to terms that contribute only at three-loop or- 
der, and with s = x in the functions on the right side. 
Combining equations (|6.42|1 . H6.43[) and (|6.44J1 gives 



d_ 
dQ' 



16tt 2 
1 



(16tt 



2 V2 



[8F 4 - 16 5 2 F 2 ]cc - 



(6.45) 



This verifies that equation (|6.4(J|) is indeed consistent 
with equations l|6.31|l and (|6.32|) . 



In the same way, one can check renormalization scale 
invariance of the complex scalar pole squared mass. One 
finds 

g * _ 4Y 2 ]z + -±^{Y 4 [8 - 20B(x, z)} 
2B(x,z) - xB(x',z) -hiz\}z (6.46) 



-16y 2 y 2 



and 



Q 



dQ 



n?2(*) 



ng(z)n^'(z)j = 20Y 4 zB(x,z) 

■2B(x,z) + xB(x',z)+Taz]z, (6.47) 

up to terms of order g 4 and terms of three-loop order and 
with s = z in all loop functions on the right side. From 
these two equations, we obtain 



-16g 2 Y 2 



Q dQ U « 



1 



■[8 ff 2 -4y 2 ]z 



1 



' (167r2)2 [ 8y4 + °(5 4 )k (6-48) 

in successful agreement with equations (|6.33() and i|6.34[) . 

Finally, for a numerical study, consider the pole 
squared masses of the neutral scalars, for a model de- 
fined by running parameters g = Y = 1 and x = 1, 
z = 0.1 at a renormalization scale Qo = 1 (in arbi- 
trary units). These choices allow the decays 4>o — > 4'+4>- 
and 4>o — > tp+'4 > - an d <fo ~~ * 0+<^--<4 an d </>o ~* ip+ip-A 
and 0q — > ip+tfi-X, and </>o — > so the pole 

squared mass will have an imaginary part correspond- 
ing to the total width. The value of equation (|6.19|l with 
(|6.21|) and (|6.22|) can be found numerically by computing 
the master integral M(z,z,z,z,x) simultaneously with 
its subordinates U(z, z, z,x), T(z,x,z), and S(x,z, z), 
[and thus V(z,z,z,x)] as described in [fj^]. The func- 
tion M(z, z, Zj_z, 0) can be computed analytically (as 
found first in [25|], and presented in the notation of the 
present paper in equation (6.27) of as can T(0, z, z), 
I(x,z,z), B(z,z), B(z',z), A{z) and A(x). I then find, 
working at Q — Qo = 1: 

\M 2 



iTM 



<po , 1-loop 



[M 2 - iTM] t 



= 0.96575- 0.03082 i, (6.49) 

, , o = 0.96423- 0.03312 i. (6.50) 

Now, to test the renormalization group scale- 
independence that was formally checked above, one can 
run the parameters x, z, Y, g from the scale Qo = 1 to any 
other scale Q using equations (|6.30|) - (|6.38|) . and then re- 
compute the complex pole squared mass in the same way. 
The results are shown in Figure The two-loop approx- 
imation to the pole squared mass features a significantly 
better scale dependence than the one-loop approxima- 
tion, as one might expect. 



VII. OUTLOOK 

In this paper, I have presented results for the two-loop 
self-energy functions for scalars in a general renormaliz- 
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FIG. 4: The real part, M 2 , and the magnitude of the imag- 
inary part, VM, of the neutral scalar pole squared mass are 
shown as a function of the logarithm of the renormalization 
scale Q, for the model described in the text. The dashed lines 
are the one-loop approximation, and the solid lines are the 
two-loop approximation. 



able theory, including the contributions of all Feynman 
diagrams that contain only one vector propagator (or 
none). This is equivalently the leading non-trivial order, 
quadratic, in gauge couplings. In some cases, notably 
Higgs scalar boson masses in the electroweak Standard 
Model and most extensions of it, these are the dominant 
two-loop contributions, because electroweak gauge cou- 
plings are smaller than the QCD and top Yukawa (and 
possibly bottom Yukawa) couplings. I have specialized 
the results obtained here to the Higgs scalar bosons of 
the Minimal Supersymmetric Standard Model, and the 
results will appear elsewhere [7i| ■ 

Further progress will require evaluating the contribu- 



tions of the remaining two-loop Feynman diagrams, con- 
taining two or more vector propagators. A hallmark of 
the strategy used here is that it is flexible, applying to 
a general renormalizable field theory. In softly broken 
supersymmctry, the full two-loop result should allow, as 
special cases, precise evaluation of the pole masses of not 
only the Higgs scalar bosons, but also the squarks and 
sleptons. 



Appendix 

This Appendix contains a few useful identities involv- 
ing the basis functions. Reference [62]| contains many 
more identities, including analytic expressions for some 
of the basis functions, and the derivatives of all of the 
basis functions with respect to s and each squared-mass 
argument, expressed algebraically in terms of the basis 
functions. 

Partial derivatives with respect to the renormalization 
scale are given by: 



Q h A[x) 



Q 8Q T{X ^ Z) 
d 

Q-QQU{x,y,z,u) 
d 

Q-QQV(x,y,z,u) 

d — 

Q-qqV(x,0 7 z,u) 



QggM(x,y,z,u,v) 



-2x, (A.l) 

2, (A.2) 

2[A(x) + A(y) + A(z) 
-x-y-z], (A. 3) 

2[A(x) + A(y) + A(z) 
-x-y-z] + s, (A. 4) 
-2A(x)/x, (A.5) 

2-2B(x,y), (A.6) 

2 + 2B(x,y), (A.7) 

-2B(x,y'), (A.8) 

{2{s + x)[B(0,x) - 1] 
+AA(x)}/(s-x) 2 , (A.9) 
0. (A.10) 



The function V(x, 0, u, z) introduced in equation (|2.2()|l 
can be rewritten explicitly in terms of the basis integrals 
as: 
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V(x, 0, z, u) = —xU (x,0,z,u)/ (s — x) 2 + 
u(z — u— x + s) 2x 



z + u 



1 



(z — u) 2 (s — x) (s — x) 2 

z(z — u + x — s) 



xT(x, z, u) 



(z — u) 3 (s — x) [z — u)(s — x) 2 



zT(z, x, u) + 



2x 



(z — u) 3 (s — x) (z — u)(s — x) 2 



uT(u, x, z) 



(z + u){S(x, z, u) - A(x)/2 +x + z + u- 5.s/8} + A(z)A(u) /[{z - u) 2 (s - x 



+B(0,x 
+ [(xz 
V(x, 0, Z, z) = 



(2x + s){A(u) - A(z)} uA(z)-zA(u) (z + u) 



A(x){A(u) - A(z)} 
(z — u)(s — x) 2 



(s — x) 2 (z~u) [z — u) 3 2(z — u) 2 

+ [(xz + xu- 2su)A(z) + (xz + xu - 2sz)A(u)] /[(z - u) 2 (s - x) 2 ], 

' ' (Az — 3a; — s)zT(z, z, x) + (s — x)xT(x, z, z) + (4z — x + s)S(x 1 z, z) — 2z/(0, z, z) 



(A.ll) 



3z(s — x) 2 

+2A(z){x -2z- A(x)} + A(x){x/2 - Az - s/2} + Az 2 - zx - x 2 + 13sx/8 + sz/2 - 5s 2 / 

-B(Q, x) [{A(z)/z + l}(s + x)/(s - x) 2 + l/6z] , (A.12) 

V(x, 0, 0, 0) = \-(s + x)T(0, 0, a;) + 25(0, 0, x) - (s + x)B(0, x) - 3A(x) +x- s/a\ /(s - xf . (A.13) 

These functions are singular as s — > x, but the functions Vffffs{ x i 0) 0j m j v ) i n ec L- l|4.44|l and Vffffv { x i 0j 0, u, v) 
in eq. (j4.83(l are smooth in that limit. 

Some identities involving vanishing squared-mass arguments are: 

U(x,y,y,0) = -T(y,0,x)+[1- A(y)/y]B(x,y) + 1, (A.14) 

U(x, 0,0,0) = -T(0,0,a:)+2B(0,a:) + l, (A.15) 

V(x, y, y, 0) = {T(0, x, y) - T(y, 0, x) - [A(y)/y + l]B(x, y) + 2[A(y) - y]B(x, y')}/2y, (A.16) 
U(x, 0, y, z) = [l/(z - y) + l/(s - x)] yT(y,x,z) + [l/(y - z) + l/(s - x)] zT(z,x,y) 

+ [2xT(x, y, z) + 2S(x, y, z) - 1(0, y, z) - A(x) - A(y) - A(z) + x + y + z- a/4] /{a - x) 

+B(Q,x)[A(y)-A(z)]/(z-y), (A.17) 
U (x, 0, y, y) = [{a -x- Ay)T(y, y, x) - AxT(x, y, y) - AS(x, y, y) + 2/(0, y, y) + 2A(x) + AA(y) 

-3x -Ay + 3s/2]/(x - a) - [1 + A(y)/y]B(0, x) (A.18) 
45(0, x, y) = (x-3y- s)T(y, x, 0) + (y - 3x - a)T(x, y, 0) + 1(0, y, y) + 1(0, x, x) - sB(x, yf 

+B(x, y)[2s + (y-x- s)A(x)/x + (x-y- s)A(y)/y] + 2[A(x) + A(y) - x - y)\ + 3s/2, (A.19) 

5(0, 0,x) = -xT(x,0,0) + (s- x)B(0,x)/2 + A(x)/2- x + 5s/8. (A.20) 



The functions in the preceding identities appear in the section Ivl quite non- unique. 

self-energy functions involving massless vector propaga- This work was supported by the National Science 
tors, making the presentation of the formulas given in Foundation under Grant No. 0140129. 
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